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Extensive density-functional calculations are performed for chemisorption of atoms in the three 
first periods (H, B, C, N, O, F, Al, Si, P, S, and CI) on the polar TiC(lll) surface. Calculations 
are also performed for O on TiC(OOl), for full 0(1 x 1) monolayer on TiC(lll), as well as for bulk 
TiC and for the clean TiC(lll) and (001) surfaces. Detailed results concerning atomic structures, 
energetics, and electronic structures are presented. For the bulk and the clean surfaces, previous 
results are confirmed. In addition, new detailed results are given on the presence of C-C bonds in 
the bulk and at the surface, as well as on the presence of a Ti-based surface resonance (TiSR) at the 
Fermi level and of C-based surface resonances (CSR's) in the lower part of the surface upper valence 
band (UVB). For the adsorption, adsorption energies -Eads and relaxed geometries are presented, 
showing great variations characterized by pyramid-shaped -Bads trends within each period. An 
extraordinarily strong chemisorption is found for the O atom, 8.8 eV/adatom. On the basis of the 
calculated electronic structures, a concerted-coupling model for the chemisorption is proposed, in 
which two different types of adatom-substrate interactions work together to provide the obtained 
strong chemisorption: (i) adatom-TiSR and (ii) adatom-CSR's. This model is used to successfully 
describe the essential features of the calculated iSads trends. The fundamental nature of this model, 
based on the Newns- Anderson model, should make it apt for general application to transition-metal 
carbides and nitrides, and for predictive purposes in technological applications, like cutting-tool 
multilayer coatings and MAX phases. 



I. INTRODUCTION 

Titanium carbide (TiC) is an intriguing material. Like 
other early transition-metal carbides (TMC) and nitrides 
(TMN) , it is characterized by a combination of properties 
that are traditionally assigned to two completely differ- 
ent types of materials: ceramics and metals. It possesses 
the physical properties (extreme hardness, high melting 
point, outstanding wear resistance, etc.) typical of ce- 
ramic materials, such as AI2O3, but exhibits at the same 
time electric and heat conductivities comparable to those 
of pure transition metals ni Such a behavior originates 
from a unique combination of covalency, ionicity, and 
metallicityi^ 

TMC's and TMN's, and in particular TiC, are of- 
ten used or suggested for use in applications where 
their interaction with the environment {e.g., growth, 
reactivity, tribology) is a crucial factor. Examples 
are: coating materials on, e.g., cemented-carbide cut- 
ting tools)^ microelectromechanical systems (MEMS),"* 
fusion-reactor wallsri biocompatible materials,® and in 
space applications jii substrate or stabilizing interlayer 
material in recording headsjS for deposited diamond 
films in electronics jiS and for the growth of carbidic 
nanostructuresii They have been suggested as good cat- 
alyst materials jiSi Also, they are one of the components in 
the MAX phases (such as Ti3SiC2 and Ti3AlC2), which 
are technologically interesting due to their unique com- 
bination of hardness and plasticity, properties that arise 
from the interaction between the layers of the MgX (for 
example, TigC) octahedra and of the A component (such 
as Si or Al) that make up their structures 

In the interaction of a material with its environment, 
chemisorption is the first step. Thus, a proper descrip- 
tion of this interaction calls for an understanding of the 



chemisorption properties and of the mechanisms behind 
them. This study focuses on TiC, since long considered 
a prototype for the early TMC's and TMN's, and on 
its (111) surface, which, although metastable, has shown 
to be the most reactive low-index surface of the early 
TMC'sii and a frequently occurring surface in the above- 
mentioned applications. 

So far, studies on the adsorption properties on the 
early TMC's have focused on the nonpolar and en- 
ergetically favored (001) surface. For the polar and 
reactive (111) surface, information is more limited. 
Studies of H and O monolayers on the (111) surface 
of TiC and of other TMC's and TMN's, as well as 
of the clean-surface electronic structures, have been 

published.l5il6,17d8,19i20,21,22,23,24,25,26,27,28 However, tO 

our knowledge, little has been done on the systemat- 
ics of the chemisorption and in particular on connecting 
these findings into a general model, in which adsorption 
of other species is taken into account. Indeed, a desir- 
able goal is to obtain a descriptive and predictive model 
for the nature of the interaction between adsorbate and 
substrate. 

For these reasons, we perform an extensive theoretical 
investigation of atomic adsorption on the TiC (111) sur- 
face, based on quantum-mechanical first-principles cal- 
culations within the density- functional theory (DFT).^— 
In order to understand the nature of the chemisorption, 
we perform a trend study, in which adsorption of H, of 
the second-period elements B, C, N, O, and F, and of 
the third-period elements Al, Si, P, S, and CI are consid- 
ered. For comparison, O adsorption on TiC(OOl) is also 
considered. In addition, in order to connect our results 
to the available experimental results, calculations on a 
full (1 X 1) monolayer of O on TiC(lll) are performed 
and the results compared to the case of atomic adsorp- 
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tion. Also, the atomic and electronic structures of bulk 
TiC and of the clean TiC(lll) and TiC(OOl) surfaces are 
studied in detail. 

In this way, the nature of chemisorption is understood 
through detailed analyses of the changes in surface elec- 
tronic structure upon adsorption and a model, describing 
the nature of the adatom-substrate interactions, is pre- 
sented. By correlating the trends in calculated adsorp- 
tion energies i?ads with the trends in calculated electronic 
structures for the different adatoms, indications about 
how £^ads depends on the adatom-TiC(lll) interaction 
are thus obtained. 

The organization of this paper is as follows. In Sec. 
II, the DFT computational method used is presented, 
including details of accuracy and convergence tests. Var- 
ious densities of states (DOS) and other key tools for 
the analysis of our calculated electronic structures are 
defined. In Sec. Ill, the results from our calculations 
are presented and analyzed. Section III. A presents the 
results for bulk TiC, with detailed emphasis on the elec- 
tronic structure. Earlier published results are confirmed. 
However, more details are given, in particular on the 
existence of direct Cp-Cp bonding states in the lower 
part of the upper valence band (UVB) of TiC. Sec- 
tion III.B describes the results for the clean TiC(lll) 
and (001) surfaces, including: surface relaxations, sur- 
face energetics, and electronic structures. Particular 
emphasis is given to the characterization of the sur- 
face electronic structures, with detailed descriptions of 
how they change compared to the bulk. Section III.C 
gives the results for the atomic adsorption on TiC (111) 
and (001), first describing in detail the calculated en- 
ergetics, atomic structures (Sec. III.C.l), and electronic 
structures (Sec. III.C. 2), and thereafter analyzing the 
obtained electronic-structure information (Sec. III.C. 3). 
The systematics of the calculated electronic structures 
are used to formulate a "concerted-coupling" model, de- 
scribing the nature of the interaction between the ad- 
sorbate atoms and the TiC(lll) surface. Section III.D 
presents the results for full (1x1) monolayer coverage of 
O on TiC(lll), including: energetics, atomic and elec- 
tronic structures, and an analysis of the applicability of 
our concerted-coupling model to the case of full mono- 
layer adsorption. In Sec. IV, we relate our electronic- 
structure results to the calculated trends in adsorption 
energies and discuss how our concerted-coupling model 
can be used in a qualitative way to understand the cal- 
culated trends in adsorption energy. In Sec. V, the main 
conclusions from the paper are briefly summarized. 



II. COMPUTATIONAL METHOD 

The calculations are performed with the DFT 
plane-wave pseudopotential code dacapo 1.30i^ The 
generalized-gradient approximation (GGA) of Perdew- 
Wang 1991 (PW91)2i is used for exchange and corre- 
lation. The atomic nuclei and core electrons are de- 



scribed with Vanderbilt ultrasoft pseudopotentials'^^ and 
the Brillouin zone (BZ) is sampled with the Monkhorst- 
Pack schemei^ Slab geometry and periodic boundary 
conditions are employed and the potential discontinuity 
in the vacuum region thus arising from the lack of slab 
mirror symmetry is corrected for.^** Full atomic relax- 
ations within the given lattice parameters are performed 
using a BFGS quasi-Newton method, as implemented by 
Bengtssoni^ On each considered system, the relaxation 
is performed until the sum of the remaining forces on 
all relaxed atoms in the unit cell is less than 0.05 eV/A, 
which corresponds to remaining forces on each individual 
atom of less than 0.01 eV/A. 

The method and the Ti, C, O, and Al pseudopotentials 
have previously been shown to yield successful results for 
alumina bulk and surface structures as well as for 
TiC/Co interface studiesi^ The pseudopotential for B 
has been succesfully employed to describe the binding in 
boron nitride. '^^ The accuracy of the other pseudopoten- 
tials used in our calculations is tested by calculating the 
atomization energies of II2, N2, F2, P2, S2, and CI2 and 
the lattice constant and binding energy of bulk Si and 
comparing these results to the experimental values and 
to the results from similar calculations. 

The electronic structures of the considered systems are 
analyzed by studying different types of local densities 
of states (LDOS), obtained from the calculated Kohn- 
Sham (KS) wavefunctions and energy eigenvalues 
En'- (i) atom- and/or orbital-projected densities of states 
[LDOS(£^)] are obtained by projecting the KS wavefunc- 
tions onto the individual atomic orbitals and plotted as 
functions of the energy relative to the Fermi level, E—Ep] 
(ii) the total density of states in real space, around a spe- 
cific energy value, [DOS(r, E)] is obtained by plotting 
|5'„(r)p, integrated over the whole BZ, for the value of 
the band index n that corresponds to the KS energy £"„ 
close to E. Also, band-structure plots are obtained by 
calculating the eigenvalues at k points located along the 
BZ symmetry lines. These plots are analyzed in detail 
by plotting |^'rik(r)P in three-dimensional real space for 
the k points located along the BZ symmetry lines. 

In addition, a measure of the localization of charge 
around individual atoms is obtained from the calculated 
DFT charge density by using the "atoms-in- molecule" - 
method of Bader.^*' The algorithm of Ref. 41 is used. In 
these calculations, the charge density is sampled on a 
grid with grid-point separation of approximately 0.07 A. 
Also, the total-electron distribution in real space is an- 
alyzed by plotting the differences between the electronic 
densities calculated for the considered systems and those 
calculated for the systems composed of the corresponding 
free atoms. 

For the calculations on bulk TiC (Sec. III. A), a plane- 
wave cutoff of 350 eV and a 6 x 6 x 6 k-point sampling 
are used. This has previously been shown to yield a good 
accuracy for bulk TiCi^ 

For the calculations on the clean TiC(lll) surface (Sec. 
III.B), a 450 eV cutoff and a 6 x 6 x 1 k-point sampling 
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are used. The TiC structure is composed of Ti and C lay- 
ers that alternate in the [111] direction (see Sec. III.B). 
In this paper, we define a "TiC bilayer" as the unit com- 
posed of one Ti and one C (111) atomic layer. Our cal- 
culations on the clean TiC (111) surface are performed on 
slabs composed of up to 15 such TiC bilayers. All consid- 
ered slabs are stoichiometric, that is, they are composed 
of equal amounts of Ti and C atoms. Therefore, one side 
of the (111) slabs is terminated by a C layer, while the 
other side is terminated by a Ti layer. Because of this, 
only the "cleavage energy" i?cieav) that is, the sum of the 
C-terminated and Ti-tcrminated surface energies, can be 
calculated: 

-Ecicav = Eccn{n) — nEb, 

where Eccwin) is the total energy of a slab with n TiC 
bilayers, as obtained from our DFT calculations, and 

E], £'ccii(n-) - £'ccii(?i - 1) 

is the bulk energy of one TiC bilayer The cleav- 

age energy corresponds to the energetical cost of cleaving 
the infinite TiC crystal into two semi-infinite crystals. 
Convergence of -Eb is obtained for n > A, disregarding 
quantum-size oscillations^^i^ which make the calculated 
cleavage energy values vary by at most ±0.3 J/m^. The 
thickness of the vacuum region is 8.75 A, corresponding 
to the thickness of four TiC bilayers. Increasing the cut- 
off energy from 350 eV to 450 eV and/or the vacuum 
thickness from four to six TiC bilayers changes the cal- 
culated £'cioav value for the unrelaxed slab by less than 
0.1 J/m^. 

For the adatom calculations on TiC(lll) (Sec. III.C), 
a slab composed of four TiC bilayers is used. In or- 
der to be able to compare our results with those from 
our study of the formation of thin films of a- and k- 
AI2O3 on TiC(lll)4^ orthogonal surface lattice parame- 
ters, resembling the structure of the K-Al2O3(001)/(001) 
surfacCjfSi are used in these calculations (see Fig.P). Each 
atomic Ti/C layer is thus composed of six atoms in a 
(2 X 3) geometry. The adatom and the top three TiC 
bilayers are allowed to relax in all directions. A 4 x 2 x 1 
k-point sampling, a 400 eV cutoff energy, and a vacuum 
thickness of at least 9.6 A are used. The adsorption en- 
ergy of the adatom is calculated as 

E^ds -^slab+adatom -^slab -^atom; 

where Ssiab+adatom IS the calculated energy of the fully 
relaxed TiC slab with the adatom, i?siab is the energy 
of the relaxed clean TiC slab, and i?atom is the energy 
of the free adatom. All values are calculated with the 
same supercell size, plane- wave cutoff energy, and k-point 
sampling. The accuracy is tested by considering the sys- 
tem with strongest adsorption energy, an O adatom ad- 
sorbed in the fee site, and repeating its relaxation with 
an 8 X 4 X 1 k-point sampling, a 700 eV plane-wave cutoff, 
a 20.6 A vacuum thickness, and a slab thickness of seven 




FIG. 1: Surface unit cell used for the adsorbate calculations 
on TiC(lll). The (4 x 3) unit cell used for the fee adatoms 
is marked by the lines and the shaded atoms. Large lighter 
circles are Ti atoms, small darker circles are C atoms, and the 
small black circle is the adsorbed atom. The present picture 
shows the relaxed structure of an adsorbed O atom in the fee 
site. The (2 x 3) unit cell, also used in the calculations (see 
text for details), is also shown (dashed arrow). 



TiC bilayers, respectively. None of these changes pro- 
duces a change in calculated adsorption energy of more 
than ±0.05 eV/adatom. 

In addition, in order to test any presence of adatom- 
adatom interactions, the relaxation of all adatoms in fee 
site is repeated using supercells in which each atomic 
layer is composed of 12 atoms in a (4 x 3) geometry (Fig. 

A 2 X 2 X 1 k-point sampling is used. Again, no 
change in adsorption energy larger than 0.05 eV/adatom 
is obtained. However, for the adatoms with large radii 
(atomic or ionic, depending on the degree of ionicity of 
the adsorbate-substrate bond, see Sec. III.C. 2) (B, F, 
Al, Si, and CI), slightly asymmetric adsorbate-substrate 
bond lengths are found after relaxation of the (2 x 3) su- 
percell. Also, evidence for adatom-adatom interactions 
are found in the calculated DOS's for these adatoms. 
Therefore, all results in this paper on the atomic adsorp- 
tion in the fee site on TiC (111) are those obtained from 
the (4 X 3) supercell. 

For the calculations on the clean and the 0-adsorbed 
TiC(OOl) surfaces (Sees. III.B and III.C), a slab com- 
posed of five atomic layers is used. Each TiC(OOl) atomic 
layer consists of equal amounts of Ti and C atoms (see 
Sec. III.B). In order to allow sufficient distance between 
the adsorbed O adatoms, the surface cell dimensions are 
chosen such that each layer is composed of 16 atoms 
(eight Ti and eight C) in a (4 x 4) geometry. A 450 eV 
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cutoff, a 4 X 4 X 1 k-point sampling, and a vacuum thick- 
ness of at least 10.9 A are used. Since the two surfaces 
of a TiC(OOl) slab are equivalent, the calculated (001) 
cleavage energy corresponds to twice the (001) surface 
energy. 

Finally, the calculations on a full 0(1 x 1) monolayer 
on TiC(lll) (Sec. III.D) are performed using both a slab 
with 15 TiC bilayers (and one atom per atomic layer) 
and a slab with 4 TiC bilayers [and six atoms per layer 
in a (2 X 3) geometry]. The same parameters are used as 
for the above-described 15-bilayer and (2 x 3) 4-bilayer 
slabs, respectively. Between these two slab, the calcu- 
lated interlayer relaxations of the O monolayer and of 
the top two TiC bilayers differ by less than 0.02 A and 
the calculated adsorption energies differ by less than 0.1 
eV/adatom. 



III. RESULTS AND ANALYSIS 
A. Bulk TiC 

1. Atomic Structure 

Bulk TiC adopts the NaCl structure. The lattice pa- 
rameter obtained from our calculation (4.332 A) IS very 
close to the experimental value (4.330 A).**^ Along the 
[111] direction, the structure is composed of alternating 
Ti and C layers, stacked in an ABC sequence (see Sec. 
III.B), with calculated interlayer distance of 1.251 A, Ti- 
C bond length of 2.166 A, and Ti~Ti and C-C distances 
within each (111) layer of 3.063 A. 



2. Electronic Structure 

The nature of the bonding in bulk TiC and 
other TMC's and TMN's has been extensively 
studied over many years, both experimentally and 
theoreticallyi^i'*'''i^^i'^^i^°i^^i^^i^^ The common viewpoint 
that emerges from these studies is that the bonding 
in several TMC's and TMN's is dominated by strong, 
directional, covalent M-X bonds (M: metal atom, X: 
C or N), which account for the high hardness of these 
compoundsi^ At the same time, a transfer of electrons 
from M to X takes place, implying a certain degree 
of ionic bonding. Further, the DOS does not vanish 
at the Fermi energy (Ep), giving the compounds a 
metallic character. Thus, the bonding in TiC is mainly 
due to covalent Ti~C bonds, in which a certain degree 
of ionicity (electron transfer from Ti to C) can be 
detected, combined with a smaller amount of metallic 
Ti-Ti bonds. Most previous studies have excluded the 
presence of direct C-C interactions a priori. How- 
ever, more recently, a COOP (crystal orbital overlap 
populations) analysis has pointed out the presence of 
bonding and antibonding C-C interactions in the lower 



and upper parts of the TiC upper valence band (UVB), 
respectively;™ 

In the following, we present our results on the elec- 
tronic structure of bulk TiC. The presence of iono- 
covalent Ti-C bonds and of metallic Ti-Ti bonds is con- 
firmed and analyzed in detail. In addition, our analysis 
clearly shows that bonding C-C states are present in the 
lower part of the UVB. Our subsequent analysis of the 
adsorption on the TiC(lll) surface (Sec. III.C.3) shows 
this to be important for the formulation of a model for 
adsorption on TiC. Therefore, we pursue in this and the 
following subsection detailed analyses of both the previ- 
ously characterized Ti-C and Ti-Ti states and the so far 
largely neglected C-C states. 

Our calculated band structure and DOS{E) for bulk 
TiC are shown in Figs. [5]and|3| The lower valence band 
(LVB), peaked at 9.5 eV below Ep, is dominated by C 2s 
states, with only small contributions from Ti. The upper 
valence band (UVB) is mainly composed of C 2p and Ti 
3d states and extends above —6.1 eV, with a major peak 
at —2.5 eV and smaller peaks at —4.0 and —1.7 eV (all en- 
ergies here and in the following are given relative to Ep). 
No C sp hybridization is detected. The main contribu- 
tion to the UVB comes from the C 2p orbitals, indicating 
the partially ionic nature of the Ti-C bond. However, the 
large overlap of Ti and C states, in particular around and 
above the main UVB peak, demonstrates a strong cova- 
lent character of the bond. The UVB and the conduction 
band (CB), on each side of Ep, form a continuum, with 
a wide, nonzero but deep, minimum in the DOS around 
Ep. The lowest value of the DOS is approximately 0.21 
states/(eV-TiC-unit), in good agreement with previous 
theoretical and experimental results4^*^ This confirms 
the metallic character of the crystal. The CB has mainly 
Ti 3d character, with a smaller C 2p contribution, which 
again indicates the partially ionic character of the Ti- 
C bond. Still, an overlap of the Ti and C states can be 
seen, which is consistent with the view that the CB arises 
mainly from the antibonding states of the Ti-C covalent 
interaction. 

Due to the octahedral symmetry around the Ti site in 
bulk TiC, the Ti 3d electrons redistribute into orbitals 
of Cg {i.e., — and z^) and t2g {i.e., xy, yz, and zx) 
symmetries The eg orbitals point toward the 
nearest C-atom sites, while the t2g orbitals point toward 
the tetrahedral holes in the fee Ti-atom network.-^^ Our 
calculated LDOS, projected onto Ti eg and t2g orbitals 
(Fig. O, shows a dominating Cg character of the UVB. 
On the other hand, both the states around Ep and in 
the CB are dominated by t2g orbitals. 

These results agree well with those obtained from 
the previous studies on the bonding nature of bulk 
Tic.2i47,48,49,50,5i,52^53 According to them, the UVB is 
dominated by bonding Ti-C a states between Ti 3d{eg) 
and C 2p orbitals {pda), with contributions from Ti-Ti 
a states between the Ti 3d{t2g) orbitals {dda) and from 
Ti-C TT states between Ti 3d{t2g) and C 2p orbitals {pdn). 

However, our results show also that (i) the Ti DOS is 
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FIG. 2: Calculated band structure (right) and total and atom-projected DOS(-E) (left) for bulk TIC. The high-symmetry points 
in the Brillouin zone are T = (0,0,0), X = ^(0, 1,0), W = ^(i, 1,0), and L = ^{\,\,\), where a is the primitive-lattice 
parameter of bulk TiC. The band-structure plot shows also which parts of the upper valence band (UVB) correspond to C-Ti, 
C-C, and Ti-Ti bonding states, respectively, according to our state-resolved analysis of the DOS (described in the text and 
illustrated in Fig. QJ. 
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FIG. 3: Calculated DOS(S) for bulk TiC, projected onto Ti d 
valence orbitals with eg and t2g symmetry, respectively, onto 
Ti s -I- p valence orbitals, and onto C s and p valence orbitals, 
respectively. 



very weak in the energetical region below the main UVB 
peak, (ii) no particular symmetry preference [s, p, d{eg), 
or d{t2g)] can be detected for the Ti states in this region, 
and (iii) C states are also present around Ep. 

In order to further investigate these issues, we examine 
real-space, three-dimensional, plots of |^„k(r)P (see Sec. 
II) for a number of states nk along all the valence bands 



of the TiC band structure. This analysis confirms the 
claims of previous studies, regarding the Ti-C and Ti-Ti 
bonds in bulk TiC. However, it also gives more detailed 
information on the types of bonding states present in 
the UVB and on where they are located in the band 
structure: (i) Ti-C bonds are mainly found along A2 
(pda bonds), A45 (pdn), Z4 (pda), Q3 and Q4 [pda and 
pd-Tr), and A34 {pdir) [Figs. EJa) and (b) illustrate this 
by the pda and pdir bonds at the minimum energy point 
of A2 and at the L34 point, respectively]; (ii) Ti-Ti dda 
bonds within the TiCjOOl} planes are found along the 
A3 band and around the point, as illustrated by the 
three-dimensional contour plot of |4'nk(r)P at the ^3 
point [Fig. He)]. 

In addition, our analysis shows the clear presence of C- 
C bonding states, mainly in the lower part of the UVB: 
the states around the X2, ^45, W23, and L2 points, and 
along the Z2, Q2, and the lower part of the A2 symmetry 
lines consist almost exclusively of C-C bonding states 
(Fig. 121). Of these states, only those along Q2 and the 
lower part of A2 contain some amount of Ti-C bonding 
states. These C-C bonding states are of both ppn (at 
A^2, ^^45, and along Z2 and the lower part of Q2 and of 
A2) and ppa (at W23, and along the higher part of Q2) 
characters. Figures ^e) and (d) illustrate this, showing 
|*nk(r)|^ at L2 and at W23, respectively. Thus, as al- 
ready suggested by the low Ti LDOS in this region, the 
lower part of the TiC UVB is almost exclusively com- 
posed of bonding C-C states, with very small Ti-C con- 






(A) (e) (f) 

FIG. 4: (Color online), (a)-(e) Three-dimensional real-space contour plots of |*I'nk(i')P in bulk TiC for the states nk at: (a) 
the point of lowest energy of the A2 band, showing Ti(eg)-C(p) a bonds; (b) the L3 point, showing Ti{t2g)-C{p) n bonds; 
(c) the X3 point, showing Ti{t2g)-Ti{t2g) a bonds; (d) the W2 point, showing C(p)-C(p) a bonds; (e) the L2 point, showing 
C(p)-C(p) TT bonds, (f) Three-dimensional contour plot of the difference between the total electron density in bulk TiC and 
the electron density of free Ti and C atoms. All subfigures are viewed perpendicular to the [111] direction, except (c), which 
shows a {001} plane. Larger (grey /lighter) balls are Ti atoms, smaller (green/darker) balls are C atoms. In subfigure (f) the 
C atoms are hidden inside the large electron clouds. Note that in subfigure (e), the Ti electron clouds lie behind and are not 
connected to the C-C electron clouds. The contour plots correspond to density values of (a)-(e) 0.17 electron states/A'' and 
(f) -1-0.02 electron states/A^. 



tributions. 

Figure^ff ) shows the difference between the total elec- 
tron density in bulk TiC and the electron density of free 
Ti and C atoms, that is, how the valence-electron dis- 
tribution is changed by the formation of bonds in bulk 
TiC. The plot confirms the partially ionic character of the 
bond, with a clear charge transfer from the Ti to the C 
atoms (this is also confirmed by an analysis of the regions 
of negative charge-density difference). At the same time, 
it can be seen that the covalent interaction is dominated 
by the Ti-C a bonds. 

Thus, we can conclude that the main contribution to 
the cohesion of bulk TiC comes from the iono-covalent 
Ti-C a bonds (in agreement with previous studies), but 
that our state-resolved analysis of the electron density 
[Figs. Ha)-(e)] shows that Ti-C tt, Ti-Ti a, C-C a, and 
C-C TT bonding states also are present in the UVB of bulk 
TiC, in different energy regions (Fig. 13): C-C bonding 
states are present in the lower part of the UVB (between 
approximately —5.8 and —3.1 eV); Ti-C bonding states 



are present from —3.9 eV all the way up to Ep; Ti-Ti 
bonding states are present between approximately —1.6 
and -1-0.7 eV. Some C-C bonding states are also found 
around the X45 point, at Ep. As will be seen later, 
knowledge about the presence and energetical location 
of these different states is important for an understand- 
ing of the mechanism behind atomic adsorption on TiC. 

Finally, we estimate the ionic character of bulk TiC 
by calculating the charge around each atom with the 
"atoms- in-molecule" method of Bader (see Sec. II). This 
yields a charge transfer of 1.5 electrons from Ti to C, in 
excellent agreement with the results from MuUiken pop- 
ulation analyses of other calculationsii^4i& 



B. Clean TiC(lll) and TiC(OOl) surfaces 

The TiC(lll) surface has been characterized, exper- 
imentally and theoretically, to be Ti(I x I) terminated 
[Figs. [SIa-b)]>i2iiSii242S Due to the partial ionicity and 
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NaCl structure of TiC, its (111) surface is polar, that is, 
it is energetically unstable due to the long-range electric 
field caused by the nonzero dipole moment, perpendicu- 
lar to the surface, of the TiC bilayer.^^ However, previous 
calculations have shown that this polarity is cancelled 
by a lowering of the ionic charge at the surface by ap- 
proximately 50%, due to the downward shift of the Ti 
antibonding CB LDOS at the surface. ^^'^^ Experiments 
have indeed shown the presence of a Ti-centered surface 
resonance (TiSR) just below Epi^ii^iSli^ 

In contrast, the TiC(OOl) surface is nonpolar, be- 
ing composed of equal amounts of positively charged Ti 
atoms and negatively charged C atoms [see Fig. El^c)]. 
Neither charge transfer nor surface resonance (SR) at Ep 
have been reported for this surface. 

In the following, we present our results from DFT 
calculations on these two surfaces. First, the relaxed 
atomic geometries and cleavage energies are presented 
and compared with previous results. Then, we focus 
on the electronic structures of the two surfaces, in par- 
ticular the (111) surface, describing in detail how they 
differ from the bulk electronic structure. These re- 
sults are in a general agreement with the few earlier 
studiesii^iiSiiLiSii^iSiiSS However, our analyses provide a 
much greater resolution of the surface electron-structure 
features and show a so far neglected wealth of differ- 
ent SR's, in particular on the (111) surface. Knowledge 
about the presence and nature of these SR's is essential 
for the subsequent analysis of the electronic structures of 
the atomic adsorbates. 



1. Geometry and Energetics 

For the (111) surface, our calculations on the 15-bilayer 
slab (see Sec. II) yield relaxations of the top six interlayer 
distances of -19.2%, +11.0%, -5.8%, +0.6%, -0.7%, 
and —1.0% [Fig. [SJb)]. The large first-layer relaxation is 
in qualitative agreement with the one observed by Aono 
et alm^ by impact-collision ion scattering spectroscopy 
(ICISS) (—30%), and all values agree well with those 
from the tight-binding calculation of Tan et alm^ The 
sum of the calculated top two interlayer distances (2.40 
A) agrees also very well with the value 2.3 A from the 
STM investigation by Ahn et alm^ The calculated relaxed 
bond length between the surface Ti and C atoms is 2.041 
A, which corresponds to a 6% decrease from the bulk 
value of 2.166 A. 

For the (001) surface, our calculations yield the rum- 
pled relaxation described in the literature, with the C 
atoms relaxing outwards and the Ti atoms inwards 
This results in a perpendicular C-Ti distance of 0.11 
A in the surface layer, in excellent agreement with re- 
cent experimental results (0.13 ± 0.04 A)m^ Our relaxed 
perpendicular distances between the surface-layer Ti and 
subsurface-layer C (2.08 A) and between C and Ti within 
the subsurface layer (0.03 A) [see Fig. Eld)] are also in 
good agreement with the experimental results (2.11+0.04 



A and 0.01 + 0.04 A, respectively).^ 

For the cleavage energies, i.e., the energy required to 
cleave the infinite crystal along a chosen plane (as de- 
fined in Sec. II), we obtain i^cioav — 11.74 J/m^ for the 
unrelaxed (111) surface and i?cicav — 11.43 J/m^ after re- 
laxation of only the Ti-terminated side of the (111) slab. 
Since the (111) slab is stoichiometric, these values corre- 
spond to the sum of the surface energies of the Ti- and 
C-terminated TiC(lll) surfaces, respectively (see Sec. 
II). However, the calculated energy gain of the relax- 
ation, 0.31 J/m^, corresponds to the surface-energy gain 
of relaxing only the TiC(lll)-Ti surface. For the (001) 
surface, a value of -Ecicav = 3.52 J/m^ is obtained after 
relaxation. Since the (001) slab is mirror symmetric with 
respect to the (001) plane, i?cicav corresponds to twice the 
surface energy for this surface. 

These values of i^cicav agree well with the results 
from previous DFT calculations performed with the VASP 
code:59 Scicav = 3.46 J/m^ for the relaxed TiC(OOl), 7.56 
J/m2 for the relaxed TiC(Oll), and 11.26 J/m^ for the re- 
laxed TiC(lll).^° Thus, we see that the energetical order 
of the three low-index surfaces of TiC is (001) < (110) < 
(111). The lower stability of (111), with respect to (001), 
is expected, since the polar nature of the (111) surface 
should decrease its stabilityi^ 



2. Electronic Structure: TiC (111) 

(a) Fermi-level surface resonance. The calculated 
DOS(i?) for the relaxed TiC(lll) surface is shown in 
Figs.Eland[7| The bulk DOS (Fig.|2Jl is recovered in the 
third TiC surface bilayer. The top two bilayers are char- 
acterized by a strong SR, peaked right beneath Ep, of 
almost exclusively Ti 3c? character (TiSR). The energet- 
ical position and Ti d character of the TiSR are in good 
agreement with the results from previous experimental 
and theoretical studies. i^iiSti^iSiiSS In addition, our or- 
bital projection of the DOS (Fig.!?)) shows that the TiSR 
can be decomposed into (i) a main peak at —0.1 eV that 
consists of mainly d(^xz,yz) and d^^y/x^-y'^) orbitals (the z 
axis is here and in the following taken perpendicular to 
the surface) and (ii) a high-energy shoulder, above Ep, 
that consists almost exclusively of d(^xz,yz) orbitals. The 
highest density of d(^xz,yz) orbitals is found at +0.3 eV. 

The real-space shape of the TiSR (Fig. ISJ reveals this 
mixed d(^xz,yz) + d(xy.x'^-y^) symmetry: the electron dis- 
tribution protrudes into the vacuum, while binding to 
neighboring surface Ti atoms. However, the TiSR elec- 
tron density is not homogeneously distributed within the 
surface plane: the TiSR avoids the regions with C atoms 
directly underneath, that is, the hep adsorption sites (see 
Sec. III.C.l). Therefore, the surface Ti-Ti bonds are 
present only in the regions corresponding to the fee ad- 
sorption sites. Also, similar plots show that, for lower 
DOS values, the TiSR density couples to states centered 
around Ti atoms lying deeper into the bulk, motivating 
our choice of calling these states SR's rather than surface 
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FIG. 5: Surface structure of TiC(lll) and (001): (a) top view of the two atomic top layers of TiC(lll); (b) the top seven 
atomic surface layers of TiC(lll), seen perpendicularly to [111], showing the ABC stacking of alternating Ti and C layers and 
the interlayer relaxations, in absolute values (d) and relative to the bulk distance (A), as calculated with our 15-bilayer slab (see 
Sec. II); (c) top view of the surface atomic layer of TiC(OOl); (d) the top four surface layers of TiC(QQl), seen perpendicularly 
to [001], showing the rumpled relaxation of this surface, as obtained from our calculations. 



states. 

This spatial extent of the TiSR agrees with the one 
suggested by Aono et alm^ and can be understood by 
considering the TiSR as containing danghng bonds that 
result from the breakage of the iono-covalent Ti-C bonds 
that cross the (111) cleavage plane. The breakage of these 
bonds causes the corresponding empty, antibonding, Ti- 
centered states of the bulk CB to collapse into more 
atomic-like orbitals, thus lowering their energy. There- 
fore, the TiSR contains unsaturated Ti orbitals that 
point toward the sites that were previously occupied by 
the C atoms, that is, the surface fee sites [Fig. IS^b)]. At 
the same time, the TiSR can also be expected to contain 
broken Ti-Ti dda bonds [Fig. EJc)] that rearrange into 
new Ti-Ti hybridizations within the surface (111) plane. 

As noted above, previous calculations have led to the 
argument that the presence of a partially-filled TiSR 
causes an electron transfer to the surface Ti atoms and a 
resulting decrease in their ionicity of approximately 50%, 
compared to the bulk,^^'^*^ Since TiC is partially ionic, 
the TiC(lll) surface is polar and should therefore be 
unstable, due to the presence of a macroscopic electric 
field caused by the nonzero perpendicular dipole moment 
of the TiC bilayer^^^ However, this field can be neutral- 
ized by inducing a surface charge that exactly cancels 
out the bulk macroscopic field. For the (111) surface 
of a crystal with the NaCl structure, such as TiC(lll), 
this surface charge should be equal to 50% of the bulk 
ionicity. ■'^^ Based on the above-mentioned calculations, it 
has been argued that the extra electronic states needed 
to accomodate this charge at the Ti-terminated surface 
are provided by the TiSRi^ 

A Bader analysis (see Sec. II) of the atomic charges 
in our TiC (111) slabs yields that (i) in the middle of 
the slabs, a charge transfer from Ti to C atoms of ap- 
proximately 1.55 electrons/atom takes place (in good 
agreement with our bulk results, see Sec. III. A. 2); (ii) 
close to the Ti-terminated (C-terminated) surface of 
the slab, there is a total net charge of approximately 



—0.8 (-1-0.8) |e|/atom. Thus, approximately 0.8 elec- 
trons/atom are transferred from the C-terminated side 
to the Ti-terminated side of the slab. This amount cor- 
responds to 0.8/1.55 = 52% of the ionicity in the middle 
of the slab, which is in excellent agreement with the pre- 
vious results of a 50% reduction in surface ionicity. 

However, a more detailed analysis of our Bader charges 
for the Ti-terminated surface shows that this surface 
charge is not located exclusively on the surface Ti atoms, 
but rather distributed among both Ti and C atoms of 
the top two TiC bilayers: the charge on the surface Ti 
atoms is +1.1 |e|/atom, that is, only ^ 0.45 more elec- 
trons than on the bulk Ti atoms. The remaining extra 
surface electrons are located on the C atoms of the first 
TiC bilayer (0.25 electrons/atom) and on the Ti and C 
atoms of the second TiC bilayer (0.03 electrons/atom on 
both Ti and C). Thus, the extra electrons do not only 
give rise to a TiSR but also an increased ionicity of the 
surface C atoms, compared to the bulk. This is char- 
acteristic of electron-structure changes also around the 
surface-bilayer C atoms, a point that is pursued further 
in the next subsection. 

Finally, it should be noted that the TiSR DOS peak is 
rather sharp, indicating a small dispersion and therefore 
a weak coupling between neighboring surface Ti atoms 
and with the bulk states at Ep. This is confirmed by our 
real-space plot of the TiSR (Fig. |SJ), which shows that 
the electron density is very weak in the surface fee sites, 
compared to the Ti sites, but also by an analysis of the 
band structure around Ep, as calculated for a slab with 
four TiC bilayers (Fig. O. Analysis of contour plots of 
|^nk(r)|^ for a number of k points along this band struc- 
ture allows an identification of the states that show a 
strong localization around the surface Ti atoms. These 
states, marked with circles in Fig. El show that the TiSR 
lies at —0.5 eV at F and disperses upwards in both the 
TM and TK directions, crossing Ep. The dispersion is 
smaller than 1 eV. Both the direction and the strength of 
the dispersion agree qualitatively with previous theoret- 
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FIG. 6: Total and atom-projected DOS(i5) for the top two 
TiC bilayers of the clean TiC(lll) surface. In the surface- 
bilayer DOS, the four UVB peaks and the peak at Ef 
are identified as C-localized (CSR1-CSR4) and Ti-localized 
(TiSR) surface resonances, as described in the text. 



ical results for TiC(lll)^^ and with experimental results 
for the largely similar NbC(lll) surface. 

(h) UVB and LVB. The UVB of the TiC(lll) surface 
bilayer extends between —5.9 and —0.7 eV and is largely 
similar to the bulk UVB. However, a detailed analysis 
reveals a number of differences: (i) the position of the 
main peak lies at —2.9 eV, that is, 0.4 eV lower than 
in the bulk UVB; (ii) the bulk-UVB high-energy peak 
(at —1.7 eV) is strongly quenched; (iii) the bulk-UVB 
low-energy peak (at —4.0 eV) is replaced by three peaks, 
at -3.7, -4.5, and -5.3 eV; (iv) the surface UVB DOS 
is more strongly localized around the C atoms than the 
bulk UVB DOS. 

Our analysis of the DOS(r, E) (see Sec. II) for all states 



within the surface UVB shows that the surface UVB can 
be divided into three regions, according to the bonding 
nature of its states: (i) a high-energy region, between 
— 1.7 and —0.7 eV, consisting of C-Ti bonding states with 
higher electron concentration around the Ti atoms; (ii) 
a middle region, between —3.3 and —1.8 eV, consisting 
of C-Ti bonding states with higher electron concentra- 
tion around the C atoms; and (iii) a low-energy region, 
between —5.9 and —3.3 eV, consisting of primarily C-C 
bonding states, with some contribution from C-Ti in- 
teractions. The high-energy region corresponds to the 
quenched high-energy bulk peak, the middle region cor- 
responds to the main peak, and the low-energy region 
corresponds to the three low-energy peaks. 

The differences between the bulk and surface UVB 
DOS's can be understood in a way analogous to our dis- 
cussion on the origin of the TiSR. The high-energy region 
of the bulk UVB is composed of bonding Ti-C states, 
with higher localization on the C atoms. Cleavage of the 
bulk structure to create the (111) surface breaks these 
Ti~C bonds and removes the corresponding C atoms. 
Therefore, these Ti-C bonding states collapse into more 
C-centered, atomic-like, states, which disappear as the 
corresponding C atoms are removed upon cleavage. This 
shows up in our surface DOS as a quenching of the bulk- 
UVB peak at -1.7 eV. 

Similarly, the formation of the surface causes the 
breakage of the C-C bonds present in the lower part of 
the bulk UVB (see Sec. III.A.2). The DOS changes in the 
low-energy UVB region can therefore be interpreted as 
evidence for the formation of C-centered dangling bonds 
or SR's (CSR's). Such an interpretation is corroborated 
by our calculated DOS(r, E) for the energetical regions 
around each of the UVB peaks in the surface DOS(i?) 
(illustrated by Fig. llO|) . A higher localization and differ- 
ent orbital symmetry of these states around the surface- 
bilayer C atoms is clearly visible. Also, at low values of 
the DOS(r, E), these states connect to the second-bilayer 
C-atom states, giving these states the status of SR's. 

Finally, the higher localization of the surface UVB 
DOS around the C atoms, compared to the bulk, is fur- 
ther evidence for the higher ionicity of the surface-bilayer 
C atoms, already discussed above in conjunction with our 
surface Bader analysis. Our analysis above shows that 
this higher ionicity is caused by the presence of CSR's on 
these atoms. 

The LVB of the TiC(lll) surface is peaked at -9.9 eV 
and is almost exclusively composed of C 2s orbitals, like 
the bulk LVB. Again, we note a shift of 0.4 eV to lower 
energy, compared to the bulk LVB. 

(c) Total electron density. The total effect of all these 
DOS changes on the surface electron distribution can be 
seen in our calculated total electron density (Fig. Illfl . 
The SR's described above are clearly visible as increases 
in electron density, compared to the bulk, above both 
surface-bilayer Ti and C atoms. 
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FIG. 7: Atom- and orbital-projected LDOS(i5)'s for the top TiC bilayer of the clean TiC(lll) surface. 




FIG. 8; (Color online). Three-dimensional real-space contour plots of DOS(r, E), summed over all occupied electronic states 
above —0.6 eV, for the clean TiC(lll) surface, showing the Fermi-level surface resonance (TiSR). The plot is viewed (a) 
perpendicularly to the surface and (b-c) from above the surface (showing only the top TiC bilayer). Plot (b) shows a higher 
DOS value than plot (c). Green balls are C atoms. The Ti atoms lie inside the electron clouds, as marked in the figures. 



3. Electronic Structure: TiC(OOl) 

For comparison, we have also calculated the DOS(i?) 
for the relaxed TiC(OOl) surface (Fig.[T21). No clear SR is 
present at Ep on this surface. The surface-layer DOS(i?) 
is very similar to the bulk DOS(i?), apart from a shift of 
the main LVB and UVB peaks by 0.8 eV toward Ep- In 
particular, the lower part of the UVB is still characterized 
by only one peak, like the bulk UVB. 

Like on TiC(lll), however, changes occur in the high- 
energy part of the bulk UVB. As discussed above, this 
part is characterized by the iono-covalent Ti-C bonds. 
Like for TiC(lll), the creation of the (001) surface breaks 
these bonds and can therefore be expected to modify 
the electronic structure at these energies. In particu- 
lar, like on TiC(lll), the strong high-energy peak of the 
bulk UVB is quenched. Also, orbital analysis shows en- 



ergetical shifts of the Tidz2 states lying in this region. 
These are, in a similar way as for TiC(lll), indications 
of electron-structure changes due to the breakage of the 
bulk Ti-C bonds that cross the (001) cleavage plane. 

Thus, compared to the (111) surface, the nonpolar 
TiC (001) surface shows no clear sign of CSR's in the 
lower part of the UVB nor of a TiSR around Ep- On the 
other hand, SR's appear in the upper part of the UVB, 
due to the breakage of Ti-C bonds. 



C. Atomic adsorption on TiC(lll) and TiC(OOl) 

To understand the nature of chemisorption on 
TiC(lll), we perform a trend study, in which adsorp- 
tion of the first-period element H, of the second-period 
elements B, C, N, O, F, and of the third-period elements 
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FIG. 9: Part of the calculated band structure around Ep for 
a TiC(lll) slab with four TiC bilayers. Real-space plots of 
|^'nk(r)|'^ for a number of the states belonging to the bands 
closest to Ef have been inspected. Those states that show 
a strong surface localization in these plots are marked in the 
figure with circles. The high-symmetry points in the two- 
dimensional Brillouin zone are F = (0,0), M = (^,0), and 
K = (|, |), where the coordinates are given relative to the 
surface reciprocal unit cell. 



Al, Si, P, S, CI are considered. For comparison, the ad- 
sorption of O on TiC(OOl) is also considered. 

For TiC(lll), all possible adsorption sites are consid- 
ered: the two different three-foldly coordinated hollow 
sites, the two-foldly coordinated bridge site, and the top 
site. The two hollow sites differ by their relation to the 
structure of the second atomic layer of the substrate: the 
site is either located directly above a second-layer C atom 
or it is not, corresponding to the hep and the fee sites, re- 
spectively. The only information present in the literature 
regarding the adsorption geometry on TiC(lll) concerns 
oxygen, which has been found to adsorb dissociatively in 
the fee site^^ and hydrogen, which has been found to ad- 
sorb in a three-fold hollow site, presumably the fee sitei2& 

For TiC(OOl), two different sites are considered: on 
top of a C atom and on top of a Ti atom. Previous 
studies have shown that oxygen adsorbs dissociatively 
on TiC(OOl), and that the preferred adsorption site is on 
top of the C atomi6ii62i63,64 

The structure of this subsection is as follows. First 
(Sec. IILC.l), the calculated energetics and relaxed ge- 
ometries of all adsorbate systems are presented. Then 
(Sec. IILC.2), the calculated electronic structures of the 
systems with adatoms adsorbed in the stable fee site of 
TiC(lll) are described and analyzed in detail. For com- 
parison, the electronic structures of the O adatom ad- 
sorbed in the hep and top sites on TiC(lll) and in the 
stable on-top-C site on TiC(OOl) are also calculated and 
analyzed. Finally (Sec. III.C.3), a model for the adsorp- 
tion mechanism on TiC(lll) is proposed, on the basis of 
these results. 



1. Energetics and geometry 

For each adatom, total-relaxation calculations are per- 
formed, yielding the ground-state equilibrium geometry 
and adsorption energy E'ads, corresponding to the en- 
ergy gained by the formation of adatom-substrate bonds 
(see Sec. II). A higher value of i^ads implies a stronger 
adatom-substrate bond. Thus, analyses of the iSads 
values give information about the preferred adsorption 
site for an adatom and a relative measure of the bond 
strengths of different adatoms on a surface. 

The calculated i?ads values for the adatoms in the dif- 
ferent sites on TiC(lll) are listed in Tabled For the fee, 
hep, and top sites, the values given are those obtained af- 
ter full relaxation of the adsorbate and substrate atoms in 
all directions, with the exception of B, C, and N adatoms 
in top site, for which small forces (~ 0.05 eV/A ) toward 
the neighboring fee sites are still present after relaxation. 
No adatom binds in the bridge site, all of them relaxing 
to the neighboring fee site. Thus, the i?ads values given 
for the bridge site are those obtained after full relaxation 
of the adatoms in only the direction perpendicular to the 
substrate and of the substrate atoms in all directions. 

All adatoms prefer to adsorb in the fee site, followed 
by the metastable hep site, where the adsorption energy 
in most cases is very close to that of the unstable bridge 
site. The bond to the top site is significantly weaker 
than to the other sites and the top-site potential-energy 
surface (PES) has a rather shallow energy mininum, a 
small deviation from the optimum relaxed position be- 
ing sufficient for destabilization and relaxation toward a 
neighboring fee site. The overall energetical preference 
for the fee site agrees with the experimental results for 
oxygen and hydrogen adsorption on TiC(lll).SSi2& 

The strongest fee bond is obtained for oxygen, almost 
9 eV/atom, followed by C, N, S, and F, the energies for 
C and N being almost the same. The weakest bonds are 
obtained for Al and H. The variation in i?ads values is 
very large, from 3.4 eV/atom for Al to 8.8 eV/atom for 
O. 

Also, the energy differences AE^ds between the differ- 
ent adsorption sites (Table vary very strongly between 
the adatoms. In particular, the bridge site can be used as 
an approximation for the transition state between the fee 
and hep sites. The energy barriers for diffusion between 
fee and hep sites can then be estimated from the differ- 
ences between the i?ads values calculated for the fcc/hcp 
sites and the bridge site. The estimated barriers for dif- 
fusion from fee to hep site thus calculated vary from 0.12 
eV for Al to 1.14 eV for C. The smallest barriers are ob- 
tained for Al, CI, Si, H, F, and S, in order of increasing 
barrier height. 

In addition, it is interesting to study the i?ads differ- 
ences between fee and hep sites. The smallest such differ- 
ences are obtained for Al, B, and Si, in order of increasing 
difference. On the other hand, the strongest preference 
for fee site is obtained for N, O, and C. 

Thus, of all considered adatoms, Al and Si appear to 
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FIG. 10: (Color online). Three-dimensional real-space contour plots of DOS(r, E) for the clean TiC(lll) surface, at selected 
energies E close to the energies of the four UVB peaks CSR1-CSR4 of the surfaces DOS(i5) (Fig.|HJ. For each UVB peak [(a)- 
(d)], two different electronic states, lying very closely in energy, are shown in order to illustrate the mixed pxy + Pz symmetry 
of each peak. The plots are viewed perpendicularly to the surface. C (Ti) atoms are green/darker (grey/lighter) balls, or lie 
inside the electron clouds, as marked in the figures. For different (same) CSR's, the plots correspond to different (same) values 
of DOS(r, E). 



TABLE I: Calculated atomic adsorption energies -Eads for the different adsorption sites on the TiC(lll) surface. For all 
adatoms, the bridge site is unstable and relaxes to the neighboring fee site. The values correspond to full relaxation of the 
adatom and surface-bilayer atomic coordinates, except for the bridge-site values, for which the adatom has only been relaxed 
in the direction perpendicular to the surface. Also given are the energy differences A_Eads between the fcc/hcp and the bridge 
sites and between the fee and the hep sites. 



-Bads (eV/atom) A^ads (eV/atom) 



atom 


fee site hep site top site bridge site 


fcc-bridf 


^e hcp-bridf 


;e fcc-hcp 


H 


3.60 


3.34 


2.31 


3.27 


0.33 


0.07 


0.26 


B 


5.68 


5.55 


3.68 


4.92 


0.76 


0.63 


0.13 


C 


7.87 


7.15 


4.69 


6.73 


1.14 


0.42 


0.72 


N 


7.86 


6.87 


4.59 


6.77 


1.09 


0.10 


0.99 





8.75 


7.93 


6.50 


7.95 


0.80 


-0.02 


0.82 


F 


6.92 


6.46 


6.01 


6.58 


0.34 


-0.12 


0.46 


Al 


3.36 


3.31 


2.78 


3.24 


0.12 


0.07 


0.05 


Si 


5.06 


4.89 


3.95 


4.74 


0.32 


0.15 


0.17 


P 


5.87 


5.61 


4.49 


5.00 


0.87 


0.61 


0.26 


S 


7.09 


6.78 


5.68 


6.75 


0.34 


0.03 


0.31 


CI 


5.48 


5.20 


4.51 


5.25 


0.23 


-0.05 


0.28 



possess the most planar PES's between the fee and hep 
sites, having both small estimated diffusion barriers and 
very similar i?ads values for the fee and hep sites. This 
indicates that the directionality of the adatom-substrate 
bond is weak for these two adatoms. This is very inter- 
esting, considering the possibility that the high plasticity 
of MAX phases such as Ti3SiC2 and Ti3AlC2 could be 
caused by a good lateral mobility between the layers of 



Si/Al atoms and TigC octahedra that make up their bulk 
structures. ^'^ 

On the other hand, the C and N adatoms show a strong 
directionality of bonds to the substrate, having both high 
diffusion barriers and a high preference for the fee site, 
with the C adatom having highest barriers. This is con- 
sistent with the strong hardness and resistance against 
deformation of bulk TiC and TiN, TiC being the hardest 
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FIG. 11; (Color online). Contour plot (through a plane 
containing the Ti and C atoms marked in the figure) of 
the difference between the total electron density at the 
TiC(lll) surface and the free-atomic densities. Green/darker 
(grey/lighter) balls are C (Ti) atoms. Only the atoms 
lying in the drawn contour-plot plane are marked by 
names, the remaining atoms lying in front of the plane. 
Red/dark (blue/black) regions correspond to positive (nega- 
tive) electron-density differences, as exemplified in the figure 
by the -|- and — signs. 



TABLE 11: Calculated, fully relaxed, atomic adsorption en- 
ergies -Eads for O in top-C and top-Ti sites on the TiC(OOl) 
surface. 

-Eads (eV/atom) 
top-C site top-Ti site 
O/TiC(001) 4.96 2.95 



one, described in the literaturei^ 

For O on TiC(OOl) (Table IHl), our results confirm the 
preference for adsorption in the on-top C site that has 
been reported in the Hterature&SSiSiSi Also, our calcu- 
lated energy difference between on-top C and on-top Ti 
site (2.01 eV/atom) compares well with the value from 
cluster DFT calculations (2.25 eV/atom).^3 

For the TiC(lll) fee site, the adsorption energies for 
second- and third-period elements show similar trends 
(Fig. 113(1 . the only difference being the very similar values 
for C and N. In both periods, the binding is strongest for 
the group- VI elements (O and S) and decreases monoton- 
ically (again, with the exception of C and N) when mov- 
ing away from group VI. Also, it decreases from second- 
to third-period elements. For O in the preferred site (on- 
top C) on TiC(OOl), the adsorption is 43% weaker than 



FIG. 12: Total and atom-projected DOS(£) for the top 
atomic layer of the clean TiC(OOl) surface. 
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FIG. 13: Calculated, fully relaxed, atomic adsorption ener- 
gies -Bads for the fee site on the TiC(lll) surface, showing 
the adsorption-strength trends. Also shown is the adsorption 
energy for O in the most stable TiC(OOl) site (on-top C). 

on TiC(lll) (Fig. 1131). 

The calculated relaxed atomic distances between the 
adsorbates in the fee site and their nearest neighbors on 
the TiC(lll) substrate (Table ^QJ show similar trends 
within the two adatom periods: both the adsorbate- 
substrate distances and the substrate Ti-Ti distances 
are smallest for group- V and/or group- VI elements and 
increase monotonically when moving away from them. 
This is again an indication of the particularly strong 
chemisorption of group- V and VI elements. For most 
adatoms, the Ti-Ti distances are smaller than the clean- 
surface values, indicating a strong adatom-Ti attraction. 
Overall, the adatom-substrate distances are larger for pe- 
riod 3, due to the extra filled electron shell in the third- 
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TABLE III: Calculated geometries of the relaxed adsorbate- 
TiC(lll) systems for all considered adatoms in fee site: 
dad-Ti and dad-c are the distances between the adatoms and 
their nearest-neighbor Ti and C atoms, respectively; ^ad-TiC 
is the perpendicular distance between the adatoms and the 
TiC(lll) surface; dTi-Ti and dc-c are the Ti-Ti and C-C 
distances within the surface TiC(lll) bilayer, close to the 
adatom; and dxi-c is the distance between the Ti and C 
atoms closest to the adatom. The bond distances in the first 
TiC bilayer on the clean TiC(lll) surface are included for 
comparison. All values are in A. 



atom 


dad-Ti 


dad-C 


^ad-TiC 


dTi-Ti 


dc-c 


dTi-C 


H 


2.01 


2.73 


1.03 


2.98 


3.08 


2.04 


B 


2.16 


3.01 


1.32 


2.95 


3.08 


2.06 


C 


1.99 


2.86 


1.10 


2.88 


3.09 


2.07 


N 


1.94 


2.80 


1.02 


2.85 


3.08 


2.06 


O 


1.98 


2.81 


1.07 


2.88 


3.08 


2.05 


F 


2.16 


2.93 


1.28 


3.00 


3.07 


2.04 


Al 


2.67 


3.54 


2.02 


3.03 


3.07 


2.05 


Si 


2.52 


3.39 


1.82 


3.00 


3.07 


2.05 


P 


2.44 


3.32 


1.72 


2.99 


3.07 


2.06 


S 


2.44 


3.30 


1.72 


3.00 


3.07 


2.05 


CI 


2.54 


3.35 


1.82 


3.06 


3.07 


2.04 


Clean surface 








3.06 


3.06 


2.04 



period adatoms. 



2. Electronic Structure 

In order to understand the nature of the chemisorption 
bond and trends, we analyze the electronic structure of 
all the considered adsorbates in the fee site on TiC(lll), 
of the O adatom in hep and top site on TiC (111), and of 
the O adatom in the on-top-C site on TiC(OOl). 

For each considered system, two different types of 
DOS's are calculated (see also Sec. II): (i) the difference 
in total, atom-, and orbital-projected DOS(£') between 
the relaxed adsorbate-f slab system and the correspond- 
ing slab with no adsorbate but with Ti and C atoms in 
the same positions [ADOS(£^)]; (ii) the distribution in 
real space of the contribution to the total DOS from spe- 
cific energy values [DOS(r, E)]. Also, the distribution of 
the total electron density around the studied adatoms, 
compared to the free-atom density, is shown. A measure 
of the localization of charge around each atom is obtained 
with the Bader method (Sec. II). 

The ADOS(i?) plots (also called adsorbate- induced 
DOS's)^^ show how the adsorption modifies the elec- 
tronic structure in the surface region, for instance, by 
the appearance of chemical bonds between the adsorbate 
and the substrate. A negative ADOS peak at a specific 
energy Eq corresponds to a quenching of the electronic 
states having energy Eq, while a positive ADOS peak 
at El corresponds to the appearance of new states at 



El. Thus, for instance, a negative peak at energy Eq, 
surrounded by two positive peaks at energies Ei < Eq 
and i?2 > Eq, can be interpreted as evidence for the hy- 
bridization of the states at Eq with adatom states, result- 
ing in the formation of bonding and antibonding states 
at El and E2, respectively. 

Projections of the ADOS(£') onto individual atoms 
and atomic orbitals in the surface-bilayer region give 
more detailed information on the nature of these bonds. 
Such information is complemented by visualizations in 
three-dimensional real space of the DOS(r, E) for all 
electronic states within the TiC(lll) UVB and TiSR re- 
gions, thus providing vivid state-resolved illustrations of 
the bonds. Finally, the Bader and total charge-density 
analyses give information on the amount of charge trans- 
fer from substrate to adatom and on the charge polar- 
ization around each adatom, thus providing indications 
on the degree of covalency vs. ionicity of the adatom- 
substrate bonds. 

Starting with the H adatom, which demonstrates 
all the key features, this subsection describes all our 
electron-structure results for the studied adatoms in deep 
detail and compares them with available experimental in- 
formation. In the following subsection (Sec. III.C.3), this 
information is analyzed and used as basis for the for- 
mulation of a model describing the adsorbate-substrate 
interaction. 

(a) H adatom in fee site on TiC(lll). The calculated 
ADOS(£') for an H adatom adsorbed in the fee site on 
TiC(lll) (Fig. ll4|) is characterized by (i) two strong neg- 
ative peaks around Ep, of almost exclusively Ti char- 
acter (one peak, of mainly d(^xy,x^-y^) symmetry, lies at 
—0.2 eV, and the other one, of almost exclusively d(^xz,yz) 
symmetry, lies at +0.4 eV); (ii) a strong positive peak 
at -1-1.1 eV, of almost exclusively Ti d character (with 
mixed d symmetry); (iii) a strong negative peak at —3.9 
eV, of mainly C character; (iv) a broad positive region 
between —6.1 and —4.1 eV (with two main peaks at —5.5 
and —5.0 eV and a smaller one at —4.3 eV), composed 
mainly of H s states, with minor contributions from the 
substrate C and Ti atoms for the two main peaks; and 
(v) a positive region between —3.6 and —0.7 eV (with 
peaks at —3.1, —2.3, and —1.4 eV), composed mainly of 
C-centered states. 

The three-dimensional, real-space, plots of DOS(r, E) 
show that (i) the peaks at —5.5 and —5.0 eV corre- 
spond mainly to strong H-C bonding states [Fig. I15r a')]. 
with weaker contributions from H-Ti bonding states; (ii) 
weaker H-Ti and H-C bonding states are present around 
the peak at —4.3 eV; (iii) bonding H-Ti states are present 
around the peaks at —3.1 [Fig. IT^ b)] and —2.3 eV. In 
addition, the DOS(r, E) summed over all occupied states 
above —0.3 eV [Fig. lTKT c')] clearly shows the depletion of 
TiSR around the adsorbed H atom. 

Our prediction of strong H s-derived peaks at —5.5 and 
—5.0 eV can be compared with the experimental finding 
of an H Is state between approximately —5.5 and —6.5 
eV for full (1 X 1) monolayer coverage of H on TiC(lll)i2i 
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FIG. 14: Calculated density of states (DOS) for an H atom 
adsorbed in the fee site on the TiC(lll) surface. The thin line 
shows the DOS projected on the H adatom. The other three 
lines show the difference in DOS ("ADOS") between the re- 
laxed slab with adsorbed fee H and a slab without H adatom 
but with Ti and C coordinates fixed to the same values as 
those of the relaxed slab+adsorbate system. The dashed (dot- 
dashed) line shows the ADOS projected on all the C (Ti) 
atoms of the first TiC bilayer of the TiC(lll) surface. The 
thick line shows the total ADOS for the adsorbate-|-first-TiC- 
bilayer system. The shaded area shows the energetical loca- 
tion of the UVB of the clean TiC(lll) surface. The dashed 
vertical lines mark the energetical locations of the four UVB 
peaks of clean TiC(lll). 

Also, our prediction of a sharp DOS decrease around Ep 
compares well with published experimental results. 

The total electron density around the H adatom (Fig. 
I16|l shows an excess of electrons (compared to a free H 
atom) around the adatom. Our Bader analysis of the 
charge localization yields an "ionicity" {i.e., charge dif- 
ference compared to the neutral atom) of — 0.64|e| (where 
|e| is the electronic charge) around the H adatom and a 
decrease of charge, compared to the clean TiC (111) sur- 
face, around its nearest-neighbor (NN) Ti and C atoms 
(Table llV)l . Thus, electrons have transferred from both 
the Ti and C atoms in the surface bilayer to the adatom. 
The extra electrons around the adatom show an almost 
spherical distribution, typical of the s orbital, although 
a polarization toward the NN Ti atom can be seen. 

(h) Second-period adatoms (B, C, N, 0, and F) in 
fee site on TiC(lll). The calculated ADOS(i;) and to- 
tal electron densities for the second-period adatoms ad- 
sorbed in the fee site on TiC (111) are shown in Figs. 1171 
and 1181 Also, the calculated ionicities (i.e., charge dif- 
ferences compared to the neutral atoms) for the adatoms 
and their Ti and C NN's, obtained from Bader analyses, 
are given in Table Hvl 

B adatom : For B, the ADOS(£;) (Fig. HT)) is charac- 
terized by (i) a strong negative double-peak structure 



TABLE IV: Charge localization around all the considered 
adatoms and their nearest-neighbor (NN) Ti and C atoms, 
as calculated from our Bader analyses. The values given in 
the table are the "ionicities" of the different atoms, that is, the 
atomic charge values compared to the neutral atoms, given in 
units of electron charge |e|. For comparison, the correspond- 
ing values for the Ti and C atoms of the surface bilayer of the 
clean TiC(lll) surface are also given. 



adatom 




ionicity 






adatom NN Ti NN C 


H 


-0.64 


+ 1.24 


-1.75 


B 


-1.09 


+ 1.20 


-1.65 


C 


-1.35 


+ 1.32 


-1.65 


N 


-1.35 


+ 1.36 


-1.68 





-1.17 


+ 1.34 


-1.68 


F 


-0.80 


+ 1.28 


-1.72 


Al 


-0.47 


+1.08 


-1.71 


Si 


-0.93 


+1.19 


-1.69 


P 


-1.12 


+1.25 


-1.68 


s 


-1.06 


+1.29 


-1.70 


CI 


-0.75 


+1.26 


-1.73 


Clean surface 




+1.09 


-1.79 



at —0.2 and +0.3 eV, of almost exclusively Ti character 
(with mainly mixed d(^xy,x^-y^) and d^^z.yz) symmetry 
for the peak at —0.2 eV, and almost exclusively d(^xz,yz) 
symmetry for the peak at +0.3 eV); (ii) a positive peak, 
of mainly "Tid(^xz,yz) character, at +1.0 eV; (iii) much 
smaller negative peaks, of mainly C character, at —4.7, 
—3.9, and —2.9 eV; (iv) a main positive peak at —0.7 eV, 
mainly composed of Ti and B p states; (v) a smaller pos- 
itive peak at —5.8 eV, of mainly B s character; (vi) much 
smaller positive peaks at —5.1 eV (of mainly B s char- 
acter) and at —3.2 eV (of mainly C and B s character); 
(vii) a positive region between —2.7 eV and the main 
positive peak, including a peak at —2.4 eV, of mainly C 
character. 

The DOS(r, E) plots show that (i) the peak at -0.7 eV 
consists of very strong B-Ti bonding states; (ii) the peak 
at —5.8 eV consists of very strong B-C bonding states 
and of weaker B-Ti bonding states; (iii) the peak at —5.1 
eV consists of weak B-Ti and B-C bonding states; (iv) 
only B-Ti bonding states are detected around —3.2 eV; 
(v) the only type of bonding states detected between —2.7 
eV and the main positive peak are weak B-Ti states; (vi) 
similarly to the case of H [Fig. llSf c)]. the large negative 
peak just below Ep corresponds to a strong depletion of 
the TiSR around the B adatom. 

The total electron-density plot [Fig. I18r a)]. together 
with the calculated Bader charges ( Table llV)l . show, sim- 
ilarly to the H adatom, a significant charge transfer from 
the surface-bilayer Ti and C atoms to the B adatom, 
yielding a B adatom ionicity of — 1.09|e|. In contrast to 
the H adatom, however, the distribution of the excess 
electrons around the B adatom shows a Pa;j,-like symme- 
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FIG. 15: (Color online). Three-dimensional contour plots of the calculated DOS(r, E), at different energies E, for an H atom 
adsorbed in the fee site on the TiC(lll) surface. All three plots are viewed from above the surface, showing only the adatom 
and the top TiC bilayer, and correspond to (a) E = —5.5 eV, showing the strong H-C bonding character of the lowest ADOS 
peak (Fig. I14II : (b) E = —3.1 eV, showing the H-Ti bonding character of the ADOS peak at —3.1 eV; and (c) the sum of all 
occupied states above E = —0.3 eV, showing the strong depletion of TiSR in the region around the H adatom. As marked in 
each plot, the smaller (green/darker) balls are C atoms, the larger (grey/lighter) balls are Ti atoms, and the red/darkest ball is 
the H adatom [in (a)-(b), the H adatom lies inside the electron cloud]. In (b), the electron cloud around the H atom connects 
only to the neighboring Ti atoms, while the C-centered electron clouds connect only to the neighboring C and Ti atoms and 
lie below the H-Ti bonding clouds. 
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FIG. 16: (Color online). Same as Fig. 1111 but for the sys- 
tem with an H adatom adsorbed in fee site on the TiC(lll) 
surface. The (darkest) red ball is the adatom. 



try, extending clearly toward the NN Ti atoms. Also, 
an electron depletion can be observed directly above the 
B adatom, again indicating a strong polarization of the 
adatom electronic distribution toward the substrate. 

C adatom : For C, the ADOS(£;) (Fig. [l7|) shows (i) 
again a large negative double-peak structure at —0.2 
and -1-0.3 eV, of predominantly Ti character (with mixed 
d(xy,x^-y^) and d(^xz,yz) Symmetry for the peak at —0.2 
eV, and only d(^xz,yz) symmetry for the peak at -fO.3 
eV), similar to the one for B but with a slightly wider 
low-energy peak; (ii) a positive peak, of predominantly 
Tid(a;z,j,z) character, at -1-1.0 eV; (iii) a main positive 



peak, of mainly adatom-C p and Ti character, at —1.3 
eV; (iv) a smaller positive peak, of mixed adatom and 
substrate character, at —3.1 eV. 

The DOS(r, E) plots show that (i) the main positive 
peak corresponds to very strong adatom-C-substrate-Ti 
bonding states; (ii) the smaller positive peak is com- 
posed of weak adatom-C-substrate-Ti and adatom-C- 
substrate-C bonding states; (iii) the large negative peak 
just below Ep is characterized by a depletion of the TiSR 
around the C adatom. 

The total electron-density plot [Fig. I18r b)] and Bader 
charges (Table llV|l are similar to those for the B adatom, 
with a significant electron transfer from the surface- 
bilayer Ti and C atoms to the C adatom, these electrons 
assuming a clear pxy-\i^e symmetry around the adatom 
and extending toward the NN Ti atoms. Also, an electron 
depletion is again found directly above the adatom, in- 
dicating a strong electronic polarization toward the sub- 
strate. Compared to the B adatom, however, the electron 
transfer is larger, the calculated C adatom ionicity being 
— 1.35|e|. Also, the polarization toward the NN Ti atoms 
is slightly smaller. 

N adatom : For N, the ADOS(£') (Fig.[l7|) displays (i) 
again a negative double-peak structure, of Ti character, 
at —0.3 and -1-0.3 eV, with the same orbital symmetries 
as for B and C, but with the low-energy peak slightly 
weaker than that for B and C; (ii) a positive peak, of 
mainly Tid^j.^ j,^) character, at -1-1.3 eV; (iii) in contrast 
to the other period-two adatoms, a broad positive region 
that extends between —5.9 and —0.9 eV, composed of five 
peaks (a main peak, at —3.4 eV, which shows mainly N p 
and Ti character, two minor peaks at —5.4 and —4.7 eV, 
which show mainly C character, and two intermediate 
peaks, at —3.9 and —2.4 eV, which show mainly N p 
character) . 

The DOS(r, E) plots show that (i) the states below 
^ —3.7 eV consist of mainly N-C bonding states, to- 
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FIG. 17: Same as Fig. I14l but for the second-period elements B, C, N, O, and F. 
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FIG. 18: (Color online). Same as Fig. 1161 but for the systems with the second-period adatoms (a) B, (b) C, (c) N, (d) O, and 
(e) F adsorbed in fee site on the TiC(lll) surface. 



gether with weaker N-Ti bonding states; (ii) the region 
around —3.4 eV consists of approximately equally strong 
N-C and N-Ti bonding states; (ii) the region above —2.9 
eV consists exclusively of N-Ti bonding states; (iii) the 
region of the negative peak just below Ep is character- 
ized by a depletion of the TiSR around the N adatom 
and by the presence of a very small N-centered dangling 
bond. 

The total electron-density plot for N [Fig.lMc)] is very 
similar to that for C. However, the adatom electrons as- 
sume a slightly more spherically symmetric distribution, 
that is, the polarization toward the substrate, and in par- 
ticular toward the NN Ti atoms, although still present, 
is smaller than in the case of the C adatom. The total 
electron transfer from the surface-bilayer Ti and C atoms 
to the N adatom is identical to that for the C adatom, 
yielding the same adatom ionicity (— 1.35|e|). 

O adatom : For O, the ADOS(£:) (Fig. [HI) shows (i) 
a negative, Ti-dominated, double-peak structure at —0.2 
and -1-0.3 eV, with the same orbital symmetries as B, 
C, and N, but with the low-energy peak weaker than 
for B, C, and N; (ii) a positive double-peak structure, 
of Ti character, above Ep (one peak, at -1-0.8 eV, shows 
mainly d(^xz,yz) symmetry, while the other peak, at -1-1.6 



eV, shows almost exclusively 6,^-2 symmetry); (iii) a main 
positive peak, of mainly O p character, at —5.5 eV; (iv) 
smaller positive peaks at —4.3 (of almost exclusively O p 
character) and at —3.5 eV (of mainly C and Opz charac- 
ter). 

The main positive peak, at —5.5 eV, is slightly broad- 
ened, with a shoulder on its high-energy side. While the 
main peak displays equally large Ti and C character, the 
higher-energy shoulder shows almost no C character. A 
projection of the peak onto O Pxy and Pz orbitals reveals 
that these are displaced from one another: the O Pxy 
peak lies at —5.4 eV, while the O pz peak lies at —5.1 
eV. The main part of the peak is thus composed of O Pxy 
orbitals, while the shoulder at higher energy is caused by 
the O Pz orbitals. 

This prediction of strong O p states around —5.5 eV 
agrees well with the experimental finding of O 2p states 
between approximately —5.0 and —6.5 eV for full (1 x 1) 
monolayer coverage of O on TiC(lll)iS& Also, the same 
experimental study indicates that the O 2pz state has 
a higher energy than the O 2pxy states, in good agree- 
ment with our results. In addition, our predicted sharp 
decrease in DOS around Ep agrees also very well with 
published experimental results>iS 
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FIG. 19: (Color online). Same as Fig. Mbi c) but for an O 
adatom adsorbed in fee site on TiC(lll). A small DOS above 
the O adatom, of clear dangling-bond character, can be seen. 

The DOS(r, E) plots for adsorbed O show that (i) the 
main, broadened, positive peak consists of very strong O- 
Ti and 0-C bonding states; (ii) only 0-Ti bonding states 
are present around the smaller positive peaks; (iii) the 
region of the negative peak just below Ep is characterized 
by a depletion of the TiSR around the O adatom and by 
a small 0-centered dangling bond (see Fig. IT^ . slightly 
more pronounced than the one for N. 

The total electron-density plot [Fig. llSf d)] shows a 
more symmetric distribution of the O adatom electrons 
than those of the previous second-period adatoms. Still, 
a polarization toward the NN Ti atoms can be seen. 
Again, the Bader charges (Table llV)l show a significant 
electron transfer from the surface-bilayer Ti and C atoms, 
yielding an O adatom ionicity of — 1.17|e|. 

F adatom : For F, the ADOS(£') (Fig. [T7|) shows (i) 
that the negative double peak around Ep is significantly 
smaller than for the previous period-two adatoms but 
still composed of a Ti (i(a;y,a:2_,y2) -|- d(^xz,yz) peak at —0.1 
eV and of a Ti d(^xz,yz) peak at -fO.3 eV; (ii) a positive 
double-peak structure above Ep, at -fO.7 and -f 1.1 eV, 
with mainly Ti d(^xz,yz) symmetry; (iii) a huge positive 
peak, of almost exclusively F p character, at —7.2 eV; 
(iv) much smaller positive peaks, of mainly C character, 
at —4.8, at —3.5, and at —3.1 eV; (v) small negative 
peaks, of C character, at —2.7 and at —2.1 eV. 

The DOS(r, E) plots show that (i) the main peak con- 
sists of mainly F-centered states, with only slight bond- 
ing to the neighboring Ti and C atoms; (ii) no F-Ti or 
F-C bonding states are found around the smaller posi- 
tive peaks; (iii) the negative-peak region just below Ep is 
again characterized by a depletion of the TiSR around the 
F adatom and by an F-centered dangling bond, stronger 
than for the previous second-period adatoms. 

The total electron-density plot for F [Fig. EJe)] is 
qualitatively different from those for the previous second- 
period adatoms. No polarization of the F adatom elec- 
tron density toward the substrate is found, rather, there 
is a polarization toward the vacuum side of the adatom. 



Still, some polarization toward the NN Ti atom can be 
detected, but this is much weaker than for the previ- 
ous second-period adatoms, the electron-density differ- 
ence being almost zero in-between the F-Ti bond. Again, 
the Bader charges (Table Hvll show a significant electron 
transfer from the surface-bilayer Ti and C atoms to the F 
adatom, yielding an almost completely ionized F adatom, 
its ionicity being — 0.80|e|. 

(c) Third-period adatoms (Al, Si, P, S, and CI) in fee 
site on TiC(lll). The calculated ADOS(£;) and total 
electron densities for the third-period adatoms adsorbed 
in the fee site on TiC(lll) (Figs. EO) and EU show strong 
similarities to those calculated for the second-period ele- 
ments. 

Al adatom : For Al, the ADOS(£;) (Fig.|2DJ) shows (i) 
a large negative double-peak structure at —0.1 eV (of 
mainly T^'id(^xy,x'^-y'^)~''~d(xz,yz) symmetry) and at -1-0.4 eV 
(of exclusively T^'^d^^z.yz) symmetry); (ii) a positive dou- 
ble peak at -1-0.9 and -1-1.6 eV, both of predominantly Ti 
character, with mixed d symmetry; (iii) a main positive 
peak, of mainly Ti character, at —0.6 eV; (iv) smaller 
negative peaks at —3.9 eV, of mainly C character, and at 
—2.8 eV, of C and Ti character; (v) small positive peaks, 
of mainly Al s character, at —5.5, —4.9, —4.3, and —3.2 
eV; (vi) a region of positive ADOS between —2.7 eV and 
the main positive peak, including a C-dominated peak 
at —2.4 eV. Orbital projection of the Al states along the 
whole UVB region reveals that they display p symmetry 
above —1.2 eV and s symmetry below —1.2 eV. 

The DOS(r, E) plots show that (i) the main peak con- 
sists of very strong Al-Ti bonding states; (ii) only Al-C 
bonding states are found in the region below —4.8 eV; 
(iii) mainly Al-Ti bonding states (with weaker Al-C con- 
tributions) are found between —4.8 and —4.2 eV; (iv) 
only weak Al-Ti bonding states are present around the 
peak at —3.2 eV; (v) almost no Al-centered states are 
found between —2.7 eV and the main positive peak; (vi) 
the negative peak just below Ep is characterized by a 
depletion of the TiSR around the Al adatom. 

The electron-density plots (Fig. I21|l show clearly the 
overall larger radii of the third-period adatoms, com- 
pared to the second-period adatoms. Despite this, simi- 
larities are present. Similarly to B, the plot for Al [Fig. 
I2ir a)] shows a strong polarization of the adatom charge 
toward the substrate and in particular toward the NN Ti 
atoms. Also, the Bader charges (Table IIVII show again 
an electron transfer from the surface-bilayer Ti and C 
atoms to the Al adatom. However, this is significantly 
weaker than for the second-period adatoms, yielding an 
Al adatom ionicity of — 0.47|e|. 

Si adatom : For Si, the ADOS(£;) (Fig. EOI) shows (i) 
a large, Ti-dominated, negative double-peak structure at 
—0.2 and -1-0.4 eV, with same orbital symmetries as for 
Al but with a slightly wider low-energy peak; (ii) a pos- 
itive, Ti-dominated, double peak at -1-0.8 eV (of mainly 
d(xz,yz) symmetry) and at -1-1.6 eV (of mixed d symme- 
try); (iii) a sharp positive peak, of mainly Ti and Si p 
character, at —1.0 eV; (iv) a smaller positive peak, of 
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FIG. 20: Same as Fig.^but for the third-period elements Al, Si, P, S, and CI. 
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FIG. 21: (Color online). Same as Fig. I16l but for the systems with the third-period adatoms (a) Al, (b) Si, (c) P, (d) S, and (e) 
CI adsorbed in fee site on the TiC(lll) surface. 



almost exclusively Si s character, at —6.6 eV; (v) much 
smaller negative peaks at —4.0 and —2.7 eV, of mainly 
C character; (vi) a small positive peak at —3.2 eV, of 
mainly C character. 

The DOS(r, E) plots show that (i) very strong Si-Ti 
bonding states are present around —1.0 eV; (ii) the states 
at —6.6 eV are almost exclusively Si centered, although 
some weak Si-Ti bonding overlap can be detected; (iii) 
almost no Si-centered states are found around the smaller 
negative and positive peaks; (iv) the negative-peak region 
just below Ep corresponds again to a depletion of the 
TiSR around the Si adatom. 

The total electron-density plot for Si [Fig. [CT bl] is 
similar to that for Al, with a strong polarization toward 
the substrate and in particular toward the NN Ti atoms. 
Compared to Al, the electron transfer from the surface- 
bilayer Ti and C atoms to the Si adatom is larger, yielding 
a Si adatom ionicity of — 0.93|e|. This is clearly visible 
in the electron-density plot, showing a higher electron 
density in the region below the Si adatom than below 
the Al adatom. 

P adatom : For P, the ADOS(i;) (Fig. OOJ) shows (i) 
again a negative, Ti-dominated, double-peak structure at 
—0.3 and +0.3 eV, with same orbital symmetries as for Al 



and Si but with weaker low-energy peak than for Al and 
Si; (ii) a positive, Ti-dominated, double peak at -1-0.8 eV 
(of mainly d(^xz,yz) symmetry) and at -1-1.5 eV (of mixed 
d symmetry); (iii) a main positive peak, slightly broader 
than for Si and of mainly P p and Ti character, at —1.6 
eV; (iv) smaller positive peaks, of mainly C character, 
at —3.1, —3.9, and —4.7 eV; (v) small negative peaks, of 
mainly C character, at —2.7 and —4.4 eV. 

The DOS(r, E) plots show (i) very strong P-Ti bond- 
ing states around the main peak; (ii) weak P-C and P-Ti 
bonding states around the peaks at —3.1 and —3.9 eV; 
(iii) weak P-Ti bonding overlap around —4.7 eV; (iv) a 
depletion of the TiSR around P and a weak P-centered 
dangling bond in the negative-peak region just below Ep. 

The total electron-density plot for P [Fig. l^ c)] is sim- 
ilar to those for the previous third-period adatoms, with 
a strong polarization toward the substrate and in partic- 
ular toward the NN Ti atoms. However, the charge den- 
sity around the P adatom is much higher than around 
Al and Si. This is reflected by the Bader charges (Ta- 
ble lIVp , yielding a stronger electron transfer from the 
surface-bilayer atoms to the P adatom than to Al and Si, 
the P adatom ionicity being — 1.12|e|. 

S adatom : For S, the ADOS(S) (Fig. OOJl resem- 
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bles that for N and shows (i) the usual negative, Ti- 
dominated, double-peak structure at —0.3 and +0.3 eV, 
with the same orbital symmetries as for Al, Si, and P, 
but with weaker peaks than for Al, Si, and P; (ii) a 
Ti-dominated double peak above Ep^ smaller than for 
previous third-period adatoms and composed of a main 
peak at -1-0.7 eV (of mainly d^^z.yz) symmetry) and of a 
secondary peak around -1-1.1 eV (of mixed d symmetry); 
(iii) a broad positive region between —5.9 and —0.9 eV, 
which, similarly to N, is composed of five peaks (at —5.4 
and —4.7 eV, which show mainly C character, and at 
—4.0, —3.3, and —2.5 eV, which show mainly S p char- 
acter) . 

The DOS(r, E) plots show that (i) the region be- 
low —3.8 eV is dominated by S-C bonding states, with 
weaker S-Ti contributions; (ii) the region around —3.3 eV 
shows equally strong amounts of S-C and S-Ti bonding 
states; (iii) the region above —3.2 eV displays only S-Ti 
bonding states; (iv) the negative-peak region just below 
Ep is characterized by a depletion of the TiSR around S 
and by a weak S-centered dangling bond, slightly stronger 
than in the case of P. 

The total electron-density plot for S [Fig. Eld)] 
shows similarities to those for the previous third-period 
adatoms. However, the electron density around the S 
adatom, although still strongly polarized toward the NN 
Ti atoms, is smaller than for the previous third-period 
adatoms. This is reflected by the smaller electron transfer 
from the surface-bilayer Ti and C atoms to the S adatom, 
compared to the previous third-period adatoms, as ob- 
tained from the Bader charges f Table HV|l . which yield a 
S adatom ionicity of — 1.06|e|. 

CI adatom : For CI, finally, the ADOS(S) (Fig. \^ 
shows (i) again a negative double peak around Ep (at 
—0.1 and -1-0.2 eV), much smaller than for the other 
third-period adatoms but with the same Ti d symme- 
tries; (ii) a positive double peak above Ep (at -1-0.8 and 
-1-1.1 eV, both with mainly Ti d(^xz.yz) symmetry); (iii) 
a main positive peak at —5.6 eV, with an almost exclu- 
sively CI p character; (iv) smaller positive peaks at —4.4 
eV (of predominantly CI p character), and at —3.5 and 
—3.3 eV (of mainly C character); (v) a small negative 
peak, of mainly C character, at —4.0 eV. 

The DOS(r, E) plots show that (i) the main peak 
consists of mainly Cl-centered states, which nevertheless 
form strong bonding coupling to neighboring C and Ti 
atoms; (ii) the states around —4.4 eV show weak Cl-Ti 
bonding coupling; (iii) all other states show little or no CI 
contribution; (iv) the negative-peak region just below Ep 
is again characterized by a depletion of the TiSR around 
the CI adatom and by a weak Cl-centered dangling bond, 
stronger than for the previous third-period adatoms. 

The total electron-density for CI [Fig. I21f e')] shows a 
significantly smaller electron density around the adatom, 
compared to the other third-period adatoms. Almost 
zero electron-density difference is found in between the 
CI adatom and its NN Ti atoms, indicating an ionic bond. 
Indeed, the Bader analysis f Table llV|) yields an almost 



completely ionized CI adatom, its ionicity being — 0.75|e|. 

(d) adatom in hep and top sites on TiC(lll). In 
order to understand whether the nature of the adsorp- 
tion is affected by the adsorption site, we also calculate 
and examine the electronic structures of the most stable 
adatom, O, in hep and top sites and compare them to 
the results for O in fee site. 

O adatom in hep site: The ADOS(£^) for an O adatom 
adsorbed in the hep site on the TiC(lll) surface [Fig. 
I22f a)] shows that (i) the negative, Ti-dominated, dou- 
ble peak around Ep (at —0.2 and -1-0.1 eV, of d^^^^y^) + 
d(xy,x^-y^) and d(xz.yz) Symmetry, respectively) is consid- 
erably smaller than in the ADOS for the fee O adatom; 
(ii) the positive, Ti-dominated, positive peak at -1-0.6 
eV (of d(^xy,x^-y^) symmetry) is also considerably smaller 
than for fee O; (iii) an additional negative peak, of mainly 
C Pz character, is present at —2.8 eV; (iv) the main 
positive peak, mainly composed of O p^y states, with 
a smaller contribution from O Pz states, lies at —4.9 eV; 
(v) smaller positive peaks are present at —4.3 and —3.7 
(of mainly Opz character) and at —2.3 and —1.6 eV (of 
mainly mixed O pz and C pz character); (vi) there is a 
shoulder on the low-energy side of the main peak, caused 
by a positive peak of C p symmetry located at —5.5 eV. 

The DOS(r, E) plots show that (i) the main peak is 
composed of strong 0-C and 0-Ti bonding states; (ii) 
the low-energy shoulder of the main peak is dominated by 
very strong 0-C bonding states; (iii) the peak at —4.3 
eV contains both 0-Ti and 0-C bonding states; (iv) 
only 0-Ti bonding states are found around the remaining 
positive ADOS peaks; (v) the region of the negative peak 
just below Ep consists of a depletion of the TiSR around 
the O adatom and of a small O-centered dangling bonds 
[Fig.Haa)]. 

O adatom in top site: The ADOS(i?) for an O adatom 
adsorbed in the top site on TiC(lll) [Fig. I22f b'l] shows 
that (i) the negative, Ti-dominated, double peak around 
Ep (at —0.2 and -f 0.4 eV, of essentially the same symme- 
tries as for hep O, although a slightly larger dz^ contribu- 
tion is present for the peak at —0.2 eV), is again consider- 
ably smaller than for the fee O adatom; (ii) the positive, 
Ti-dominated, peak at -1-0.9 eV (of mainly d(^xz,yz) sym- 
metry) has a similar height as the corresponding peak 
for fee O; (iii) the main positive peak (of mainly O Pxy 
character) lies at a considerably higher energy than for 
fee and hep O, at —3.3 eV; (iv) other positive peaks are 
present at —5.5 and at —4.7 (of mainly C character), 
and at —2.5 cV (of mainly O Pz character); (v) there is a 
shoulder on the low-energy side of the main peak, caused 
by a C-centered peak located at —4.0 eV. 

The DOS(r, E) plots show almost exclusively 0-Ti 
bonding states. Still, a clear 0-C coupling is present 
at the low-energy shoulder of the main peak, at —4.0 
eV, and a weaker 0-C coupling is present at the main 
peak, at -3.3 eV. The sum of the DOS(r, E) for all 
occupied states above —0.3 eV [Fig. EST b)] shows that 
the top O adatom causes a depletion of the TiSR in the 
three neighboring fee sites. As opposed to fee and hep O, 
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FIG. 22: Same as Fig. I14l but for an O atom adsorbed in (a) the hep site and in (b) the top site on the TiC(lll) surface. 





FIG. 23: Same as Fig. llSf c'l but for an O atom adsorbed in (a) the hep site and in (b) the top site on the TiC(lll) surface. 
The contour plots correspond to the same value of the DOS(r, E) as the contour plot in Fig. 1191 



however, the TiSR is non- vanishing directly below the O 
adatom, due to the presence of the nearest-neighbor Ti 
atom. 

(e) adatom on TiC(OOl). For comparison, the elec- 
tronic structure for O on TiC(OOl) is also examined. 

The ADOS(i?) for an O adatom adsorbed in the on- 
top-C site of the TiC(OOl) surface (Fig. |^ is character- 
ized by (i) a strong negative peak, of predominantly C 
character, at —1.5 eV; (ii) a positive double-peak struc- 
ture, of mainly O and C character, at —6.3 and —5.9 
eV (orbital projection reveals that the more bound peak, 
at —6.3 eV, consists of O and C states with exclusively 
Pxy symmetry, while the peak at —5.9 eV shows only pz 
symmetry); (iii) smaller positive peaks, of mixed adatom 
and substrate character, at —3.7, —2.8, and —2.2; (iv) an 
0-dominated positive peak at —0.6 eV; (v) small positive 



peaks, of mainly Ti character, above Ep (at -1-0.2, -fO.7, 
and +1.4 eV). 



3. Analysis of the Electronic Structure: The 
Concerted- Coupling Model 

In this subsection, the above-described electron- 
structure results for the atomic adsorption on TiC(lll) 
are analyzed and interpreted in terms of the Newns- 
Anderson (NA) chemisorption model. This yields a 
model for the nature of the chemisorption on TiC(lll) 
that we call the "concerted-coupling model" , since it de- 
scribes the adatom-substrate interaction as arising from 
the concerted action of two contributions: (i) a strong, 
in the NA sense, coupling between the adatom state (s) 
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FIG. 24: Same as Fig. I14l but for an O atom adsorbed in the 
on-top-C site on the TiC(OOl) surface. 

and the substrate TiSR, and (ii) a coupling, whose 
strength varies between the different adatoms, between 
the adatom state(s) and the substrate CSR's. 

First [Sec. III. C. 3(a)], the essential points of the NA 
model are briefly summarized, providing a theoretical 
framework for our subsequent analysis of the electron- 
structure results for fee H on TiC(lll) [Sec. III.C.3(b)]. 
A model for the nature of the H-TiC(lll) interaction 
is thus proposed. Then [Sees. III.C.3(c)-(f)], this model 
is succesfully applied to the description of our electron- 
structure results for the other systems described in Sec. 
III.C.2. Thus, we conclude that this model is suited for 
a general description of the nature of chemisorption on 
TiC(IlI). In Sec. III.D, the model is also shown to yield 
a good description of the chemisorption of a full (1x1) 
monolayer of O on TiC(lll). In Sec. IV, the energeti- 
cal consequences of the model are discussed in a quali- 
tative way. By correlating the trends in electronic struc- 
tures with the trends in calculated adsorption energies, 
we point out which features of the concerted-coupling 
model are important for the chemisorption strength on 
TiC(lll). 

(a) The Newns- Anderson (NA) model for chemisorp- 
tion. In the NA model, chemisorption on metal surfaces 
is described with a Hamiltonian originally introduced 
for the description of magnetic impuritiesS The solu- 
tions to the NA Hamiltionian provide a theoretical frame- 
work for the description of the basic mechanisms behind 
chemisorption on metal surfaces. It has been used, for ex- 
ample, to successfully describe chemisorption trends on 
transition- metal surfaces (the d-band model) i^^iZSiii 

Here, the basic concepts arising from the NA model 
are presented. These are then used in our analysis of the 
electron-structure results for chemisorption on TiC(lll). 

Traditionally, two limiting cases are considered in the 
solution of the NA Hamiltonian: "weak" and "strong" 



chemisorption. ^^'^^ 

In the weak limit, the coupling matrix element between 
adatom and substrate states is small compared to the 
width of the substrate band. A single solution is thus 
obtained, E « e*-|-A(£*), where e* is the effective atomic 
level of the adatom (that is, the free atomic level shifted 
due to the electrostatic interaction of the adatom with 
the metal surface), and A(e*) is small compared to the 
substrate band width. If E falls within the energetical re- 
gion of the substrate band, the resulting adsorbate level 
is broadened. Otherwise, the resulting DOS consists of a 
single state outside the substrate energy band, together 
with a typically small and continuous contribution ex- 
tending over the whole substrate band. 

In the strong limit, the adatom-substrate coupling ma- 
trix element is large compared to the substrate band 
width. This yields three solutions to the NA Hamilto- 
nian: two single states, one on each side of the substrate 
band, together with a weak continuous part extending 
over the whole substrate band. These two single states 
can be regarded as the "bonding" and "antibonding" 
states, respectively, of a "surface molecule" formed by 
the adatom and its neighboring substrate atoms. 

By applying the NA model in a tight-binding approxi- 
mation, Ref. 68 shows that these bonding and antibond- 
ing states appear gradually with increasing size of the 
ratio P' 1 13, where (3' is the adatom-substrate hopping in- 
tegral and (3 is the substrate hopping integral. At small 
/?'/ (3 values, the coupling is weak and the adsorbate level 
is broadened. As (3' / (3 increases, this level broadens fur- 
ther and splits gradually into two separate levels. 

In addition to the weak and strong chemisorption, it is 
therefore reasonable to speak also of an "intermediate" 
chemisorption, in which the adsorbate level is strongly 
broadened and split into bonding and antibonding sub- 
peaks. The stronger the chemisorption is, the more are 
these subpeaks separated from one another. 

(h) H adatom in fee site on TiC(lll). The ADOS(£;) 
for fee H on TiC(lll) fFig. I14|) shows a strong negative 
peak in the energetical region of the TiSR. This indi- 
cates that the adsorption of H causes a strong quenching 
of the TiSR around the adatom, a fact confirmed by our 
real-space visualization of the states in this energetical 
region [Fig. ITCT c')]. Thus, the TiSR participates in the 
formation of adatom-substrate bonds. Since the TiSR is 
rather sharp, having weak dispersion and being well lo- 
calized around the surface Ti atoms [Sec. III. B. 2(a)], the 
coupling between the adatom state and the TiSR can 
be expected to be strong, in the NA sense, and should 
therefore give rise to well separated bonding and anti- 
bonding surface-molecule states (Fig. [23) • The bonding 
state should lie below the atomic level of the adatom and 
the antibonding state should lie above the TiSR energy. 

Indeed, the ADOS(i?) shows the appearance of a pos- 
itive peak above the TiSR and of a broad positive re- 
gion below it. The peak above the TiSR can thus be in- 
terpreted as the antibonding adatom-TiSR state, while 
the positive region below the TiSR can be interpreted as 
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FIG. 25: Schematic diagram of the concerted coupling between the adatom state and the TiC(lll) Ti- and C-centered surface 
resonances (TiSR and CSR1-CSR4). 



bonding adatom-TiSR states. 

At the same time, the adatom level can also be ex- 
pected to hybridize with the substrate UVB states, in 
a way similar to the weak coupling to the sp band 
on transition-metal surfaces On TiC(lll), however, 
such an interaction should be stronger than on pure 
transition metals, due to the more localized nature of 
the TiC(lll) UVB, compared to the transition-metal sp 
band. Indeed, our ADOS(i?) shows a strong negative 
peak at the energetical location of one of the TiC(lll) 
UVB peaks, indicating a strong coupling between this 
peak and the adatom state. As shown in Sec. III. B. 2(b), 
this UVB peak (called CSR3 in Fig. 0) corresponds to 
surface-localized C states, or CSR's. Indeed, the negative 
ADOS peak is mainly localized around the surface-layer 
C atoms. 

Furthermore, a closer inspection of the ADOS reveals 
that negative, C-centered, peaks are also present at the 
energetical locations of the two lowest TiC(lll) UVB 
peaks, CSRl and CSR2, and that the energies of the 
positive peaks lie in-between the energies of the UVB 
peaks. Thus, this indicates that the adatom level couples 
rather strongly with the TiC(lll) CSR's, giving rise to 
bonding and antibonding states that combine to form the 
observed positive ADOS peaks fFig.l^. 

According to this interpretation, the low-energy ADOS 
peaks should be characterized by bonding adatom-C 
states, while the high-energy peaks should lack such 
states. Our DOS(r, E) plots confirm this: at lower ener- 
gies, bonding adatom-C states dominate over adatom-Ti 
states [Fig. ll5r a')]: as the energy increases, this adatom-C 
character diminishes in favor of the bonding adatom-Ti 



states; above the main negative ADOS(i?) peak, only 
adatom-Ti states are found [Fig. lTsT b)]. 

Our electron-structure results for fee H on TiC(lll) 
indicate therefore that the chemisorption arises from 
the concerted action of two different types of adatom- 
substrate couplings: (i) a strong coupling between the 
adatom state and the TiSR and (ii) a somewhat weaker 
coupling between the adatom state and the CSR's that 
are present in the low-energy part of the TiC(lll) UVB. 
We call this a concerted coupling. This mechanism has 
similarities to the d-band model for chemisorption on 
transition metals. ^^•^''■^^ Here, however, the weak cou- 
pling to the sp band is replaced by a stronger coupling to 
the more localized TiC(lll) CSR's, increasing the impor- 
tance of this coupling for the chemisorption strength. In 
particular, the following analyses show that the strength 
of this coupling varies significantly between the different 
adatoms and can therefore be expected to play a role in 
the chemisorption-strength trends. 

In order to analyze these variations and trends, we 
now turn to the analysis of the electronic structures for 
the second- and third-period adatoms and of whether the 
concerted-coupling model is valid also for these adatoms, 
which are characterized by frontier orbitals of p symme- 
try, in contrast to the simpler s symmetry of the H or- 
bital. 

(c) Second-period adatoms in fee site on TiC(lll). 
Like for H, the ADOS's for the second-period adatoms 
in fee site on TiC(lll) (Fig. [T7|) show a strong nega- 
tive peak in the energetical region of the TiSR, together 
with positive peak(s) just above the TiSR and regions of 
sharp or broadened positive peaks below it. Thus, also 
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for the second-period adatoms, there is a strong couphng 
between the adatom state and the TiSR. The area of the 
negative ADOS peak increases shghtly as B ^ C and 
decreases asC— >N— >0— >F, indicating corresponding 
increase and decreases, respectively, of the adatom-TiSR 
couphng strength. 

The energetical region above the TiSR looks generally 
similar for all second-period adatoms, with a Ti-centered 
positive single- or double-peak around ~ -1-1.0 eV. Like 
for H, this can be interpreted as the antibonding solution 
of the adatom-TiSR coupling. On the other hand, the 
positive region below the TiSR, which can correspond- 
ingly be interpreted as composed of bonding adatom- 
TiSR states, differs between the different adatoms. For 
all adatoms, it contains a positive region characterized 
by states of mixed substrate-Ti d, substrate-C p, and 
adatom p symmetries. This region is successively shifted 
to lower energies as the adatom nuclear charge Z in- 
creases. This shift agrees with the interpretation of this 
region as the bonding solution of the adatom-TiSR cou- 
pling. Such a bonding solution must lie below the energy 
of the free adatom state. As Z increases, the valence 
electrons of the free adatom are more strongly bound to 
its nucleus, that is, their energies are lowered. As the 
free-adatom state decreases in energy, so must also the 
bonding solution of the adatom-TiSR coupling. 

As Z increases, the bonding adatom-TiSR state over- 
laps therefore with different parts of the substrate elec- 
tronic structure, at successively decreasing energies. In 
particular, our calculated ADOS's show that: for B, it 
lies at the upper edge of the substrate UVB region and 
is very sharp; for C, it lies in the high-energy region of 
the UVB and is slightly broadened; for N, it lies in the 
middle of the UVB region and is strongly broadened and 
split into subpeaks; for O, it lies in the low-energy region 
of the UVB and is broadened and somewhat split into 
subpeaks; and for F, it lies well below the UVB region 
and is very sharp. In addition, our DOS(r, E) plots show 
that, for all adatoms, only adatom-Ti bonding states are 
present in the high-energy region of the UVB, while both 
adatom-Ti and adatom-C bonding states are present in 
the middle- and low-energy regions of the UVB, and that 
the adatom-C contribution increases successively as the 
energy decreases. 

Thus, the strength and nature of the coupling between 
the adatom and the TiC(lll) UVB states vary signifi- 
cantly between the different adatoms. For B, the sharp- 
ness of the positive ADOS peak indicates that it only cor- 
responds to bonding adatom-TiSR states [as confirmed 
by the DOS(r, E) plots] and that no coupling takes place 
with the UVB states. 

For C, the positive ADOS peak is slightly broadened, 
indicating a weak coupling of the adatom state to the 
states in the upper region of the UVB. Since these states 
are Ti-dominated [see Sec. III. B. 2(b)], no strong adatom- 
C bonds are formed. However, a small subpeak is formed 
just below the main UVB peak, containing very weak 
adatom-C bonding states. This indicates a very weak hy- 



bridization with the C-centered states of the main UVB 
peak. 

For N, the strong broadening and partial splitting 
of the positive ADOS peak indicate a rather strong 
adatom-UVB coupling. Like for the fee H adatom, 
positive peaks are formed in between the peaks of the 
substrate- UVB DOS. At the same time, negative C- 
localized peaks are present at the energies of the sub- 
strate UVB peaks, indicating that the CSR's in the sub- 
strate UVB hybridize with the adatom state. The for- 
mation of subpeaks in the positive ADOS can there- 
fore be interpreted as the formation of bonds between 
the adatom and the C atoms in the surface TiC bilayer. 
Similarly to the case of the H adatom, these subpeaks 
correspond to combinations of bonding and antibonding 
adatom-C states. This is confirmed by the DOS(r, E) 
plots, which show a successively stronger dominance of 
bonding adatom-C states with decreasing energy, as well 
as the absence of any bonding adatom-C states at ener- 
gies higher than the main substrate-UVB peak (CSR4). 

At the same time, the differences between the N and 
H ADOS's can be attributed to the different symme- 
tries of the adatom valence states (s for H and p for N). 
For example, the much stronger coupling to the TiSR 
for N (as well as for the other second- and third-period 
adatoms) than for H (illustrated by the much stronger 
negative ADOS peaks around Ep for N) can be un- 
derstood from the better spatial overlap of the TiSR 
dxz,yz + dj.y x^_y2 orbitals with thep (mainly, the p^y) or- 
bitals of the second- and third-period adatoms than with 
the s orbital of the H adatom. 

For O, the splitting of the positive ADOS peak indi- 
cates again an adatom-UVB coupling. Like for N, posi- 
tive subpeaks are present in between the substrate-UVB 
peaks, while negative, C-localized, peaks are present at 
the energies of the UVB peaks. As described in Sec. 
III. C. 2(b), the main, broad, O p peak, which straddles 
the lowest-energy UVB peak, is actually composed of an 
O Pxy peak and of an O Pz peak, each lying on opposite 
sides of the lowest-energy UVB peak. Therefore, also for 
O, the ADOS can be interpreted as evidence for rather 
strong hybridizations between the adatom state and the 
CSR's. In particular, the main positive ADOS peaks lie 
here in the UVB region that is dominated by the CSR's, 
that is, below the main peak (CSR4) of the substrate 
UVB. As our DOS(r, E) plots show, it is in this region 
that the strongest bonding adatom-C states are found. 

It is also interesting to notice that the most strongly 
bonding 0-C states have O Pxy symmetry. This agrees 
with the experimental findings for full (1x1) monolayer 
coverage of O on TiC(lll) [see Sec. III.C.2(b)].28 Appar- 
ently, the stronger bonding of the O Pxy orbitals, com- 
pared to the O Pz orbitals, is due to the geometry of the 
fee site, which favors lateral (xy) interaction with both 
the Ti and the C atoms of the TiC(lll) surface bilayer. 

Finally, for F, the sharpness of the positive ADOS peak 
shows that no adatom-UVB coupling occurs, due to the 
fact that the adatom level lies well outside the energetical 
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region of the UVB. 

Thus, hke for the fee H adatom, two different types 
of couphngs can take place also for the second-period 
adatoms: one between the adatom state and the sub- 
strate TiSR and one between the adatom state and the 
substrate CSR's. However, their strengths vary signifi- 
cantly with varying Z within the adatom period. In par- 
ticular, the strength of the coupling to the CSR's varies 
from no coupling (B and F) to rather strong coupling (O 
and N). 

Finally, we note that, for B, the B s state lies in the 
low-energy region of the substrate UVB. This causes a 
slight broadening of the B s level, as well as the formation 
of a weak subpeak slightly above the main B s peak. The 
DOS(r, plots show that, like for the adatom -p states, 
this is due to a hybridization between the adatom s state 
and the CSR's of the surface TiC bilayer. 

(d) Third-period adatoms in fee site on TiC(lll). 
The calculated electronic structures for the third-period 
adatoms in fee site on TiC(lll) show trends that are very 
similar to those for the second-period adatoms. Over- 
all, however, two main differences can be detected: (i) 
the negative ADOS peak is smaller than for the second- 
period adsorbates and (ii) the positive ADOS peak be- 
low the TiSR lies at higher energies than for the second- 
period adsorbates. Both differences are due to core or- 
thogonalization, (i) resulting from the larger radii of the 
third-period adatoms, which cause smaller overlaps be- 
tween adatom and TiSR wavefunctions and thus weaker 
adatom-TiSR couplings and (ii) from the higher energy, 
within a group, of the valence electron state of a free 
third-period atom. 

Still, however, the two types of couplings described 
above are present. Like for the second-period adatoms, 
the adatom-TiSR coupling strength appears to slightly 
increase as Al Si and to decrease as Si — > P — > S 

CI. The coupling between the adatom p state and the 
CSR's follows also a similar trend as for the second-period 
adatoms, however, with differences due to the higher free- 
adatom energy levels. The coupling to the UVB is now 
very weak for both Al and Si, both with positive ADOS 
peaks around the top of the substrate-UVB region. The 
ADOS for P resembles that for C, the slightly lower en- 
ergy of the P peak, however, causing a slightly stronger 
adatom-UVB coupling, shown by the slightly stronger 
peak broadening and the formation of a low-energy sub- 
peak with slightly stronger adatom-C bonding character. 
In the same way, the ADOS for S resembles that for N, 
with strong subpeaks formed by the hybridization be- 
tween the adatom state and the CSR's. 

For CI, the main positive ADOS peak lies in the low- 
energy region of the substrate UVB, at a slightly lower 
energy than that of the O main peak. Thus, the CI state 
might be expected to behave similarly to the O state 
and couple strongly to the CSR's. However, our calcu- 
lated ADOS's show a smaller broadening and subpeak 
formation for CI than for O. Further, the DOS(r, E) 
plots show that the CI peak is composed of states that 



are strongly localized around the CI adatom, with a very 
weak coupling to substrate atoms. The adatom-UVB 
coupling is thus much weaker for CI than for O. The 
reason is clearly seen in our Bader calculations of the 
atomic charges f Table HVj) . For all adatoms, electrons 
are transferred from the surface-bilayer Ti and C atoms 
to the adatom. For the group- VII elements (F and CI) 
this transfer is large enough to almost completely fill the 
outer valence-electron shell of the adatom, thus almost 
completely ionizing these adatoms. This is confirmed by 
our analysis of the total electronic density around each 
adatom (see Sec. III.C.2 and Figs. El and I21|l . which 
shows an increase in electron density along the adatom- 
Ti bond for all adatoms except for F and CI. This is 
indicative of covalent adatom-Ti bonds for all adatoms 
except for F and CI, which show a prevalently ionic type 
of bond. The almost filled valence-electron shell of CI is 
therefore almost inert toward further hybridization with 
the CSR's, resulting in a weaker adatom-UVB coupling 
than in the case of O. 

Thus, also for the third-period adatoms, chemisorp- 
tion is characterized by the concerted action of the two 
different types of coupling: adatom-TiSR and adatom- 
CSR's. Again, within the period, the coupling strengths 
vary significantly with Z . In addition, they differ be- 
tween the two periods, due to the extra electron shell in 
period three. 

Finally, it can be noted that for Al, the adatom s state 
lies inside the energetical region of the substrate UVB 
and is therefore broadened and split into subpeaks of 
adatom-C and adatom-Ti bonding character. 

(e) O adatom in hep and top sites on TiC(lll). In 
order to test whether the concerted-coupling model can 
also be used to understand our results for adsorption in 
other sites than fee, we now analyze the electronic struc- 
tures for the O adatom in hep and top sites. 

For both hep and top O adatom, the ADOS's show 
negative peaks in the energetical region of the TiSR, to- 
gether with positive peaks just above and broader pos- 
itive regions below. Thus, the coupling to the TiSR is 
active in these adsorption sites as well. However, for 
both hep and top O, the resulting bonding states {i.e., 
the positive ADOS region below the TiSR) lie at higher 
energies than for fee O , implying a weaker adatom-TiSR 
coupling. 

For the hep adatom, the reason for this can be found in 
the symmetry of the TiSR. As shown in Sec. III. B. 2(a), 
due to its dangling-bond character, the TiSR extends to- 
ward the fee site and avoids the hep site. Thus, the over- 
lap between the adatom states and the TiSR is smaller 
in the hep site than in the fee site, resulting in a weaker 
hopping integral /?' between the adatom and its NN Ti 
atoms. 

This is confirmed by our DOS(r, E) plots for the en- 
ergetical region of the TiSR, which show that also for 
the hep O adsorbate [Fig. I^ a)] , the TiSR depletion oc- 
curs around the fee site close to the adatom. Also, the 
relaxation displaces the hep adatom slightly toward the 
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neighboring fee site, indieating that the adatom-TiSR 
hybridization is stronger in the fee site. 

As described in Sec. III. C. 3(a), the degree of broaden- 
ing of the adsorbate level and of its splitting into bond- 
ing and antibonding states depends on the strength ra- 
tio /3'//3, where /? is the substrate hopping integral.^'* A 
smaller value of (3' / (3 implies a smaller splitting. Thus, 
the weaker adatom-TiSR overlap in the hep site, com- 
pared to the fee site, causes a smaller downward shift of 
the resulting bonding-state energy from the free adatom- 
state energy. Indeed, the main positive peak below the 
TiSR in the hep ADOS hes 0.6 eV higher than in the fee 
ADOS. Because of this, the coupling with the UVB states 
is also changed, the hep O state coupling more strongly 
to the CSR2 peak of the substrate UVB, compared to 
the fee O state. 

In the top site, the adatom-TiSR overlap is larger than 
in the hep site, implying a higher value for /?'. How- 
ever, the atomic coordination is lower in the top site. 
The smaller splitting of the bonding and antibonding 
adatom-TiSR states can therefore be understood from 
the smaller value of the second moment of the adatom 
LDOS, /i^^^ — NP', where N is the coordination of the 
adatom. 

The difference between the top and fee ADOS's is 
much more pronounced than that between hep and fee, 
the main positive peak of the top ADOS lying 2.2 eV 
higher than in the fee ADOS. Again, the interaction 
with the substrate UVB is considerably changed because 
of this, the strongest coupling now being to the main 
UVB peak (CSR4). Correspondingly, the adatom-UVB 
coupling is more localized around the surface-bilayer Ti 
atoms than the C atoms. This is also due to the site 
geometry, which makes overlap with C-localized states 
negligible. Still, though, some clear 0-C coupling is de- 
tected in the lower part of the main positive ADOS peak, 
showing that the adatom-CSR coupling is weakly active 
also for top O. 

We conclude that the concerted-coupling model is ca- 
pable of explaining the electronic-structure features for 
adsorption in hep and top sites as well. Like for the fee 
site, the chemisorption in hep and top sites can be de- 
scribed as composed of two interacting components: an 
adatom-TiSR coupling and an adatom-CSR coupling. 

(f) O adatom on TiC(OOl). Since the nonpolar 
TiC(OOl) surface lacks the TiSR characteristic of the po- 
lar (111) surface [Sec. III.B.3], a comparison between the 
electronic structures of O adsorbed on the two surfaces 
can provide some information on the importance of the 
TiSR for the adsorption. 

Our calculated ADOS for O in the most favorable 
TiC(OOl) site (on-top C) shows the presence of only one 
negative peak fFig. I24|l . which lies just above the main 
peak of the substrate UVB and is mainly localized around 
the substrate C atom. Thus, only one type of coupHng is 
present here, between the adatom state and C-localized 
states in the upper part of the substrate UVB. As dis- 
cussed above [Sec. III.B.3], this is the only UVB region 
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FIG. 26: Structure of TiC(lll)-0(l x 1): (a) top view of the 
three top surface layers; (b) the top six surface layers, seen 
perpendicularly to [111], showing the stacking sequence and 
the interlayer relaxations, in absolute values (d) and relative 
to the TiC bulk distance (A), calculated with our 15-bilayer 
slab (see Sec. II). 

that shows signs, in our results, of SR's. The ADOS 
shows that this coupling forms bonding peaks just below 
the UVB region, of mixed O and C character, and an 
antibonding region around and above Ep. 

Thus, the weaker adsorption on TiC(OOl), compared 
to TiC(lll), can be understood to arise from a lack of 
coupling to Ti-centered SR's on (001). 

D. Oxygen (1 x 1) monolayer on TiC(lll) 

So far, experimental studies of adsorption on TiC(lll) 
have primarily been concerned with full oxygen (1 x 1) 
monolayer cover age sSSiSS Therefore, we conclude our in- 
vestigation with calculations on the O monolayer on 
TiC(lll). As for the case of atomic adsorption, first, 
the adsorption energetics and geometry are investigated. 
Then, the electronic structure is analyzed in order to un- 
derstand the bonding nature. In particular, we consider 
to what extent our concerted-coupling model for atomic 
chemisorption can be used in the situation of full mono- 
layer coverage. 

1. Energetics and geometry 

On TiC(lll), full monolayer coverage of 0(1 x 1) 
has been observed experimentally J-^ The unreconstructed 
0(1 X 1) structure is reported to be stable up to ap- 
proximately 1000°C, where it starts to reconstruct to 
(V3x V^)-R30°. The 0(1 x 1) atoms are found to occupy 
the hollow sites with no C directly underneath, that is, 
the fee sites [see Fig. ECT a)].^''^ 

Our calculated adsorption energies for 0(1x1) in both 
hollow sites (fee and hep) and in top site f Table FV)! agree 
with the experimental results of the fee site being the 
stable one. 

According to our calculations, the TiC(lll) surface re- 
laxation is significantly affected by the fee O overlayer. 
The four top interlayer relaxations of TiC(lll) (as cal- 
culated with our 15-bilayer slab, see Sec. II) are now: 



29 



TABLE V: Adsorption energies iJads for a full monolayer of O 
(1 X 1) on TiC(lll), as calculated with our 15-bilayer slab (see 
Sec. II). Adsorption in fee, hep, and top sites is considered. 

-Bads [eV/atom] 
fee site hep site top site 
8.12 7.31 5.48 



+5.1%, +0.2%, -0.2%, and +0.9% [Fig. E^b)], as com- 
pared to the values -19.2%, +11.0%, -5.8%, and +0.6% 
for the clean TiC(lll) surface [Fig. (Sfb)]. The distance 
between the O layer and the TiC(lll) surface is 0.891 
A. The 0-Ti bond is thus quite strong, pulling up the 
top Ti layer away from the second C layer and toward the 
chemisorbed O layer. Also, we note that the relaxation is 
limited to the first surface layer (the interlayer distances 
change by less than 0.02 A from the second TiC bilayer 
on). 

The ICISS study of Souda et al?^ reports about the 
presence of two different O sites, with different Ti-0 in- 
terlayer distances (1.0 A and 0.8 A, respectively). No 
sign of this is present in our results, our calculated in- 
terlayer distance lying in the middle of these two values. 
We stress here that the convergence of our calculations 
has been carefully tested and that full relaxation using 
both a (1 X 1) and a (2 x 3) surface unit cell has been 
performed [see Sec. II for details]. In addition, the calcu- 
lation with the (2 x 3) cell has been performed both with 
and without symmetry constraints on the calculated KS 
wavef unct ions . 

Our relaxed 0-Ti and 0-C bond lengths are 1.98 A 
and 2.83 A, respectively. 



2. Electronic Structure 

Like for the case of atomic adsorption, we now turn 
to the analysis of the electronic structure for full (1x1) 
monolayer of O on TiC (111), in order to understand the 
nature of the bond. 

In particular, the calculated ADOS for fee 0(1 x 1) 
(Fig. I77|l is compared with that for the fee O adatom 
(Fig. [T7jl . For fee 0(1 x 1), we note that (i) the negative 
peak around Ep is broadened and consists of a single 
peak, at +0.1 eV, of almost exclusively Ti d(^xz,yz) and 
d{xy,x^-y^) character; (ii) the positive peak in the low- 
energy part of the UVB is broadened and split into two 
subpeaks, at —5.2 and —4.4 eV, of almost exclusively 
O character (the smaller peak, at —5.2 eV, has almost 
exclusively O Pxy symmetry) ; (iii) smaller positive peaks 
are present at —3.4 eV (of mainly O Pz symmetry), at 
—2.1 eV (of mainly O Pz and C Pz symmetry), and at 
— 1.1 eV (of mainly C Pz symmetry). 

The DOS(r, E) plots for the fee O (1 x 1) monolayer 
show that (i) the peak at —5.2 eV is composed of O- 
O, 0-Ti, and 0-C bonding states; (ii) the peak at —4.4 
eV is composed of 0-Ti and 0-C bonding states; (iii) 
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FIG. 27: Same as Fig. El but for a full (1 x 1) monolayer of 
O atoms adsorbed in the fee sites of the TiC(lll) surface. 

the remaining positive peaks contain only 0-Ti bonding 
states. 

These results agree with experimental studies that 
show that the O Pxy orbitals lie at lower energy than 
the O Pz orbitalsiS^ 

3. Analysis of the Electronic Structure 

Our calculated ADOS (Fig. shows a strong deple- 
tion of states in the energetical region of the TiSR. This 
indicates that the adatom-TiSR coupling is active also 
in the case of full monolayer (1x1) adsorption of O on 
TiC(lll). In addition, the DOS(r, E) plots show the 
presence of adatom-C bonding states in the low-energy 
region of the UVB, implying that the adatom-CSR cou- 
pling is also still present. Thus, our concerted-coupling 
model appears to be valid also for the case of full mono- 
layer coverage of O on TiC(lll). 

However, the calculated ADOS for 0(1 x 1) adsorption 
shows a considerably stronger broadening of the adatom 
p peak than for atomic O fFig. [T7|) . as well as a splitting 
of the adatom p peak into two subpeaks. The DOS(r, E) 
plots show that 0-0 bonding states are present in the 
lower subpeak and absent in the higher subpeak. Thus, 
the broadening and splitting of the p peak arise from 
the coupling between neighboring O atoms within the O 
monolayer, which gives rise to the formation of bonding 
and antibonding 0-0 states. The bonding 0-0 states 
dominate in the lower part of the O p peak, while anti- 
bonding 0-0 states dominate in the upper part of the O 
p peak. 

It is interesting to notice that the O Pxy states are more 
important for the most bonding states, while the O Pz 
states are more important at higher energies. Obviously, 
the O Pxy states overlap more strongly with the neighbor- 
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ing O atoms in the monolayer, thus providing stronger 
binding. However, the O Pxy states are more bonding 
also in the case of atomic adsorption (as discussed in 
Sees. III.C.2 and III.C.3) and are therefore more favor- 
able also for the coupling with the TiC(lll) surface. This 
can be understood by considering that the O adatom lies 
in the fee site, that is, in the middle of the two triangles 
formed by the NN Ti and C atoms. Thus, the interaction 
with the NN Ti and C atoms is predominantly lateral, as 
exemplified by Figs. Iisr al-fbl. and as discussed also in 
Sec. III.C.3(c). 



IV. DISCUSSION 

In the previous sections, a concerted-coupling model is 
presented to describe the mechanisms behind chemisorp- 
tion on the polar TiC(lll) surface. The model is based 
on detailed analyses, within the theoretical framework 
of the Newns- Anderson (NA) model, of the electronic 
structures for the clean surface, for first-, second-, and 
third-period adatoms in fee site, for the O adatom in hep 
and top sites, and for a full (1 x 1) monolayer coverage of 
O. Also, for comparison, O adsorption on the nonpolar 
TiC(OOl) surface is considered. 

According to this model, chemisorption on TiC(lll) 
is caused by the combination of two different types of 
adatom-substrate couphngs: (i) a strong, in the NA 
sense, coupling between the adatom state(s) and the Ti- 
localized SR (TiSR) of TiC(lll) and (ii) a coupling, 
whose strength varies significantly between the different 
adatoms, between the adatom state(s) and the TiC(lll) 
UVB states, in particular the C-localized SR's (CSR's) 
in the low-energy part of the UVB. 

At the same time, our calculations yield characteris- 
tic, pyramid-like, trends for the atomic fee adsorption 
energies -Eads within each adatom period, with strongest 
adsorption for the group- VI elements (O and S). In par- 
ticular, the adsorption energy for fee O adatom is very 
high, almost 9 eV/adatom, reflecting the high reactivity 
of the polar TiC(lll) surface. 

In this section, we discuss how the concerted-coupling 
model can provide an understanding of these adsorption- 
energy trends and of the extremely large variations of ad- 
sorption energies between adatoms (from 3.4 eV/adatom 
for fee Al to 8.8 eV/adatom for fee O) and between sur- 
face orientations [the i?ads value for O on TiC(OOl) being 
43% smaller than on TiC(lll)]. 

In particular, we wish to address the following ques- 
tions: (i) How can we understand the -Eads trends within 
each adatom period? (ii) How can we understand the 
lower i?ads values for third-period adatoms? (iii) Why 
is chemisorption strongest in the fee site? (iv) Why is 
chemisorption stronger on TiC(lll) than on TiC(OOl)? 

In the following, the trends in our calculated electronic 
structures for the considered systems are qualitatively 
correlated to the calculated -Bads trends. In this way, 
we identify features in the concerted-coupling model that 



can be used to understand the calculated -Bads trends. 

In order to answer question (i), we start by exam- 
ining the second-period adatoms. As discussed in Sec. 
III. C. 3(c), as the adatom nuclear charge Z increases 
within the period, the two components of the concerted 
coupling vary, in different ways. 

The variations in area of the negative ADOS peaks 
around Ep indicate variations in the adatom-TiSR cou- 
pling strength: slightly increasing as B — > C and strongly 
decreasing as C ^ N — > O — > F. The strong decrease from 
C to F can be understood as a consequence of the decreas- 
ing energetical overlap between the adlevel and TiSR en- 
ergies, which decreases the strength of the adatom-TiSR 
coupling. The slight increase from B to C is based on a 
widening of the negative ADOS peak due to the down- 
ward shift of the main positive peak. Thus, it may simply 
be due to a cancellation effect between the positive and 
negative peaks of B, which partially overlap. 

Thus, the trend in chemisorption strength that arises 
from the adatom-TiSR coupling can be qualitatively il- 
lustrated as in Fig. r28f a') , where we allow an uncertainty 
in the B ^ C trend. 

On the other hand, the total strength of the adatom- 
CSR couplings increases as the adlevel energy approaches 
the energy of the CSR's, that is, the region between the 
main substrate-UVB peak (CSR4) and the lower edge of 
the substrate UVB. The best energetical overlap between 
the adlevel and the CSR's is obtained for the O adatom, 
for which all positive ADOS peaks of adatom p symmetry 
below Ep lie in the CSR region. 

Thus, the trend in chemisorption strength that arises 
from the adatom-CSR coupling can be qualitatively il- 
lustrated as in Fig. [2HIb). 

As can ben seen, the adatom-TiSR trend docs not 
agree with the Bads trend, while the adatom-CSR docs. 
We conclude that, in order to reproduce the calcu- 
lated trend in adsorption energy for the second-period 
adatoms, a description in which both the adatom-TiSR 
and adatom-CSR couplings are present is needed. In 
particular, our analysis points toward the novel result 
that the C-localized states in the lower part of the TiC 
UVB, although less important for the bonding nature in 
bulk TiC, play an important role in the description of the 
nature of chemisorption on the TiC(lll) surface. 

We also note that the almost constant value of Bads 
between C and N can be explained by the fact that the 
two coupling-strength trends go in opposite directions as 
C ^ N. 

For the third-period adatoms, the variations in 
coupling-strength trends follow similar trends to those 
described above for the second period [see Sec. 
III. C. 3(d)]. The main difference is due to the higher en- 
ergy of the third-period free-adatom level, which causes 
a better overlap with the CSR's for CI than for S, in con- 
trast to our calculated B^ads trend. As discussed in Sec. 
III. C. 3(d), however, a weakening of the Cl-CSR couphng 
is introduced by the fact that the CI adatom is almost 
fully ionized. 
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FIG. 28: Qualitative graph of the suggested trends in adsorption energy Bads for the second-period adatoms caused by (a) the 
adatoiri-TiSR coupUng and (b) the adatom-CSR couphngs. In (a), the uncertainty in the B ^ C trend (see text for details) 
is marked by an error bar. 



Finally, wc comment on the presence of couplings with 
the CSR's for the B and Al s states [described in Sees. 
III.C.3(c)-(d)]. In the same way as for the coupling be- 
tween the CSR's and the p states of CI, it should be 
expected that a coupling between CSR's and adatom s 
states is weak, due to the filled nature of the s electron 
shells for B and Al. 

Question (ii) can now be simply understood. As dis- 
cussed in Sec. III. C. 3(d), two main differences, both 
due to core orthogonalization, can be detected between 
the electronic structures of the second- and third-period 
adatoms: overall smaller negative ADOS peaks around 
Ep for the third-period adatoms and overall higher ener- 
gies for the positive ADOS peaks below Ep for the third- 
period adatoms. The first difference indicates a weaker 
adatom-TiSR coupling for the third-period adatoms. 
The second difference causes, within an adatom group, 
a weaker energetical overlap between adlevcl and CSR's 
for the third-period adatoms. Thus, both the coupling to 
the TiSR and to the CSR's are weaker for period-three 
adatoms, in agreement with our calculated E^ds trends. 

Question (iii) can also be answered through the 
concerted-coupling model. As discussed in Sec. 
III. C. 3(e), the preference for fee site over hep site is due 
to the spatial extension of the TiSR, which extends to- 
ward the fee site and avoids the hep site. The result- 
ing smaller spatial overlap between hep adatom state 
and TiSR, compared to the fee adatom, causes a weaker 
adatom-TiSR coupling. This in turn causes a smaller 
downward shift of the bonding adatom TiSR state and 
consequently also a weaker coupling between the adatom 
state and the CSR's. Thus, both the adatom-TiSR and 
adatom-CSR couplings are weaker in the hep site than 
in the fee site. 



For the top adatom, the lower adatom coordination 
causes a considerably smaller downward shift of the 
adatom-TiSR bonding state than for both hep and fee 
adatoms, again wc^akening the chemisorption strength. 
Also, due to the top-site geometry, the spatial overlap 
between the top-adatom state and the C-loealized states 
is negligible, thus strongly weakening also the adatom- 
CSR coupling. 

Finally, question (iv) can be simply understood [Sec. 
III.C.3(f)] from the lack of a TiSR on the TiC(OOl) sur- 
face. Thus, one of the couplings present on the TiC(lll) 
surface is here absent, causing a considerably weaker 
chemisorption. 

In summary, the concerted-coupling model can be used 
to successfully explain the essential features of the cal- 
culated i^ads trends. In particular, the coupling of the 
adatom state to the C-localized states (CSR's) of the 
TiC(lll) surface is expected to play an important role 
for the chemisorption trends on this surface. This is qual- 
itatively different from the case of, e.g., adsorption on 
transition- metal surfaces, where the weak coupling to the 
sp band does not affect the adsorption trends. 



V. CONCLUSIONS 

Chemisorption of atoms in the first three periods in 
the table of elements (H, B, C, N, O, F, Al, Si, P, S, and 
CI) on the polar TiC(lll) surface, as well as of atomic 
O on the nonpolar TiC(OOl) surface and of full (1x1) 
monolayer coverage of O on TiC(lll) are studied in a 
systematic way by density-functional theory. Adsorption 
energies and geometries arc calculated in fully relaxed 
configurations. Detailed analyses of the electronic struc- 
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tures of the different systems, as well as of bulk TiC and 
of the clean TiC (111) and (001) surfaces reveal the nature 
of chemisorption on this scientifically and technologically 
interesting material. At the same time, new detailed in- 
formation is obtained on the electronic structures of bulk 
and surface TiC, showing, in particular, the presence of 
direct C-C interactions and of C-localized surface reso- 
nances and their importance for the adsorption mecha- 
nisms. 

The adsorption on TiC(lll) is found to be more com- 
plex and versatile than on pure metal surfaces. The 
results for the atomic adsorption energies i?ads of the 
second- and third-period elements show strong varia- 
tions, characterized by pyramid-shaped trends within 
each period. An extraordinarily high binding energy is 
obtained for the O atom, 8.8 eV/adatom, followed by C, 
N, S, and F. The lowest energies are obtained for Al, H, 
and Si. 

To describe the variety in calculated adsorption ener- 
gies, we propose a concerted-coupling model that involves 
the concerted action of two different types of couplings: 

(i) between the adatom state(s) and the Ti-centered 
Fermi-level surface resonance (TiSR) of TiC(lll) and 

(ii) between the adatom state(s) and the C-centered sur- 
face resonances (CSR's) lying in the lower part of the 
TiC(lll) upper valence band (UVB). 

This model explains successfully the essential features 
of the calculated trends in i?ads- For example, the ex- 
ceptionally strong chemisorption of O is understood to 
arise from this combination of two different adsorption 
mechanisms. 

The strong variations in adsorption strengths show the 
versatility of TiC for technological applications. The 
strong O adsorption makes TiC(lll) a good substrate for 
oxide-layer growth, a fact that is used in multilayer CVD 
(chemical-vapor deposition) alumina/Ti(C,N) coatings 
on wear-resistant cutting tools. At the same time, the 
high Eads value for S suggests it to be a strong candi- 
date in the competition for adsorption sites, which can 
explain the presence of S found in the CVD alumina coat- 



ings under certain conditions, affecting the alumina phase 
composition^^ 

Although proper bulk and growth calculations obvi- 
ously are needed, our calculated diffusion-barrier esti- 
mates might be indicative for the synthesis and prop- 
erties of the MAX phases. Their good plasticity relates 
to a weak bonding between the (111) face of TigC layers 
and, e.g., Al or Si.^'^ Our diffusion-barrier estimates for 
Al and Si (of the order of 0.1-0.3 eV, corresponding to 
activation temperatures of ~ 50-130 K)^^ suggest good 
lateral mobility between the Al/Si and TigC layers. 

The fundamental nature of the concerted-coupling 
model, based on the Newns- Anderson model, should 
make it apt for general application to transition-metal 
carbides and nitrides. In separate publications, the ad- 
sorption features and electronic structures of TiN(lll) 
are compared to those for TiC(lll), showing the applica- 
bility of the concerted-coupling model also for adsorption 
on TiN(lll)ii^*Z& For instance, higher diffusion-barrier 
estimates are found on TiN(lll), which might provide 
ideas for the fine-tuning of the structural composition of 
the MAX phases. 

Further studies addressing these questions would be 
highly interesting. There are many questions that still 
remain to be answered. 
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Note that this TiSR has been called surface state in the 
literature, but we prefer here to adopt the proper term 
surface resonance (SR) , since its energy lies in the pseudo- 
gap between the bulk UVB and CB, that is, in a region of 
nonzero bulk DOS. 



